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Uber das elektrostatische Feld in Metallen, 
insbesondere an der Oberflache gegen einen Nichtleiter 


Von 
FRANZ OLLENDORFF 


(Mitteilung aus dem Forschungsinstitut fiir Elektrotechnik des Technion, 
: Israel Institute of Technology, Haifa) 


Mit 12 Textabbildungen 
(Eingegangen am 28. Marz 1960) 


Inhaltsiibersicht: 1. Ziel der Arbeit — 2. Grundlagen — 3. Beweglichkeit und Diffusion der freien Elek- 
tronen — 4. Das elektrostatische Feld im eben begrenzten Halbkérper — 5. Numerische Angaben — 6. Uber 
die Natur der Austrittsarbeit — 7. Das Kathoden-Grenzfeld — 8. Der ideale Kondensator — 9. Zusammen- 
fassung. 


1. Ziel der Arbeit 


In der klassischen, FARADAY-MAXWELILschen Elektrodynamik wird das elektrische Verhalten 
jedes homogenen, isotropen Stoffes durch die gleichzeitigen Angaben seiner [relativen] skalaren 
Dielektrizitatskonstanten ¢ und seiner ebenfalls skalaren elektrischen Leitfahigkeit x erschépfend 
beschrieben. Insbesondere ist der Vektor E der jeweils herrschenden elektrischen Feldstarke 
mit dem ebendort auftretenden Vektor 7 der elektrischen Leitungsstromdichte durch das differen- 
tielle Ohmsche Gesetz 

ee E (1) 
linear verkniipft. Durch das Symbol runder Klammern das skalare Produkt zweier Vektoren 
bezeichnend, finden wir aus (1) die spezifische Leistung 

p = (j E) = x(EP, (2) 
welche laut Aussage des JOULEschen Gesetzes je Raumeinheit des durchstrémten Stoffes und je 
Einheit der laufenden Zeit in Warme umgesetzt wird. 

Auf den stationaren Zustand des Systemes spezialisierend, kénnen wir auf Grund der dann vom 
Induktionsgesetz verbiirgten Wirbelfreiheit des elektrischen Feldvektors E diesen gewi8 als 
negativen Gradienten eines elektrischen Skalarpotentiales darstellen 


E =—gradp. (3) 


Verscharfen wir jetzt die bisherige Annahme des nur stationaren Zustandes zur Voraus- 
setzung eines elektrostatischen Feldes, so zieht dessen konservativer Charakter die Forderung 


p = x(E)* = 0 (4) 
nach sich; sie enthalt folgende Alternative: 
a) Ein Feld endlicher Starke kann nur in Isolatoren auftreten 
LO ME ==10: 5 (5) 
b) Jede selbst noch so geringfiigige elektrische Leitfahigkeit fihrt zum Zusammenbruch 


der elektrostatischen Feldstarke 
ee tO Tete 0), (6) 


so daB gema8 (3) im Innern jedes Leiters ein ortsunabhangiges Potential herrscht: 
y = const. (7) 
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9 _ F. Otrenvorrr: Uber das elektrostatische Feld in Metallen eae 


Lassen sich diese wohlbekannten Schliisse mit unseren modernen Auffassungen tiber die Struktur 
der Materie in Einklang bringen? 

Beziiglich der Isolatoren lehrt die Erfahrung die Existenz von Stoffen, deren Eigenschaften 
der allgemeinen Beschreibung (5) sehr nahe kommen; insbesondere reprasentiert der leere 
Raum einen allerdings nur ideellen ,,Stoff‘‘ dieser Art, mit welchem wir nichtsdestoweniger 
weiterhin bedenkenlos operieren diirfen und wollen. 

Dagegen sind uns unter den vielfachen Leitern der Elektrizitaét keine Stoffe bekannt, welche 
sich in annehmbarer Genauigkeit entsprechend den einfachen Voraussagen (6) und (7) verhalten, 
und von jeher bilden die hierher gehérigen Fragen der sogenannten Kontaktpotentiale ein 
Schmerzenskind der Elektrostatik. Uberdies wiirde die supponierte Existenz der Gleichungen (6) 
und (7) an der Grenze eines Leiters gegen einen Isolator schwerwiegenden Einwanden ausgesetzt 
sein: Ungeachtet der dort von der Wirbelfreiheit des elektrischen Feldvektors £ verlangten Stetig- 
keit des Potentiales y wiirde ja dort vermége (5) die normal zur Grenzflache weisende Kompo- 
nente £, der elektrischen Feldstarke um einen endlichen Betrag springen. Eine solche Unstetig- 
keit kann zwar mittels des Begriffes einer flachenhaft ausgebreiteten Ladung unschwer quantita- 
tiv erfaBt werden, doch ist diese Konzeption physikalisch unhaltbar, da sie dem lediglich mathe- 
matisch ausfiihrbaren GrenzprozeB zu einer unme8bar hohen Raumladungsdichte entstammt. 

,,90 schwer es uns auch werden mag, uns von alten, langvertrauten Vorstellungen geistig zu trennen: 
Unter dem Zwange der vorgenannten Erkenntnisse miissen wir uns dazu verstehen, die Gleichungen (6) und (7) 


als wissenschaftliche Grundlage einer modernen Elektrostatik der Leiter zu verwerfen. Auch im Leiter ist das 
elektrostatische Potential y in der Regel ortsabhangig, so daB dort eine endliche elektrische Feldstarke herrscht.“ 


In der vorliegenden Arbeit stellen wir uns die Aufgabe, diese gewiB einschneidende Behaup- 
tung an dem vielleicht einfachsten Beispiel der Leiter zu beweisen: An den Metallen, indem wir 
deren als unzureichend gebrandmarkte, klassische Beschreibung durch eine in sich logisch ein- 
wandfreie ersetzen, welche an Hand passender Modellvorstellungen zumindest einem Teil der 
Beobachtungstatsachen gerecht wird. Uber die prinzipielle, wissenschaftliche und didaktische 
Bedeutung dieser Frage hinaus darf ihre Behandlung eine durchaus zeitgemaBes technisches 
Interesse beanspruchen: Die Antwort erschlieBt uns den Zugang zu einer sachgemaBen Unter- 
suchung von Vorgangen der thermischen und der Feldemission, sie fiihrt uns zum Verstandnis der 
Wirkungsweise elektrisch regelbarer Kondensatoren, Kristallgleichrichter und ahnlicher Gerate 
der Festkérper-Elektronik, sofern man nur erforderlichen Falles die hier auf Metalle begrenzten 
Uberlegungen sinngema8 auf Halbleiter erweitert. 


2. Grundlagen 
Wir beschaftigen uns mit dem elektrostatischen Feld innerhalb eines chemisch reinen Metalles 
vom Atomgewicht A und dem spezifischen Gewicht y, welchem wir makroskopisch die Gestalt 
eines festen, rechtkantigen Blockes der Lange a, der Breite 6 und der Hohe c zuschreiben. Un- 
geachtet seiner jeweils tatsachlich vorliegenden Mikrostruktur werden wir diesen Quader fortan 
mit einem idealen, st6rungsfreien Einkristall identifizieren; in ihm bilden die Kerne der Metall- 
atome eine regelmaBiges Raumgitter, innerhalb dessen sich die Metallelektronen entsprechend 
ihren BLocuschen Funktionen mehr oder minder frei bewegen kénnen. Um uns von der von Fall 
zu Fall vorgegebenen Sonderform der Gitterzellen zu befreien, gehen wir von dem wahren Metall 
zu einem Modell iiber, welches sich durch ein kubisches Gitter auszeichnet. Mit dem Symbol 


V0 602 = 1078 (8) 
die AVOGADROsche (LoscumipTsche) Zahl der Atome je Kilogrammatom bezeichnend, miBt dann 
ees 
A 
d= ae (9) 


den Abstand je zweier Nachbaratome jenes Modellgitters oder, kurz gesagt, seine Gitter- 
konstante. An Hand dieses Begriffes verscharfen wir die Angabe der Kantenlangen a, 6 und c 
zur Voraussetzung ganzzahliger Verhaltnisse 


a b Cc 
ages Me ara Caemahe (10) 
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so daB die Grenzflachen des Metallblockes mit je einer Gitterebene des Modellkristalles koinzi- 
dieren. 

Innerhalb des Metalles einschlieBlich seiner Umgebung orientieren wir uns an Hand eines ruhen- 
den Bezugssystemes der KARTESIschen Koordinaten x, y und z, welche beziehentlich parallel zu 
den Kanten a, 6 und c gerichtet seien; sein Ursprung O mége zunachst entsprechend Bild 1 im 
Zentrum des Kristalles liegen. Sei nun 


Ae tee ot, ey Ss (11) 
so kOnnen wir innerhalb des Metallblocks zwei Zonen wesentlich verschiedener physikalischer 
Eigenschaften unterscheiden: 


a) Das Zentralgebiet werde durch die Angaben 
1<<|5/<5a; 1 <2 )<in; 1K Zi<i» (12) 
Va) 


definiert. In ihm wird der zeitliche Erwartungswert des elektrischen Feldes — und nur von 
diesem kann im Rahmen der verlangten, statischen Auffassung des Systemzustandes die Rede 
sein — durch ein beziiglich x, y und z dreifach periodisches, p 

zeitfreies elektrisches Skalarpotential g je der einheitlichen, 
primitiven Wellenlange d beschrieben, welches als Mikro- 
potential bezeichnet werden moge; wir kOnnen es sonach 
analytisch bereits durch die abzaéhlbare Gesamtheit seiner 
dreifach unendlich vielen Fourierkoeffizienten erschépfend 
darstellen, welche dem jeweils vorliegenden Metalle — nach 
seiner ,,Stilisierenden‘‘ Anpassung an unser Modell — von der 
Natur seiner Atome eingepragt sind und daher dieses ein- fiir 
allemal charakterisieren. 


x 
d 


Da das Mikropotential in die Schroedingergleichung der 
Elektronenwellen eingeht, regelt es deren Kinetik mittels 
der Gesamtheit der physikalisch -zulassigen Lésungen, 
welche, in hier zureichendem Umfange, durch die Stich- Didi Onienticomeia Facial 
worte ,,BLocHfunktionen“ und ,,Energiebander“ gekenn- 
zeichnet werden mégen. Beiendlichen Werten derim Innern des Metalles herrschenden absoluten 
Temperatur wird die Periodizitat des Mikropotentiales durch die dann geweckten Schwingungen 
des Raumgitters stochastisch gestért, und die hieraus resultierenden Modulationen der im Metall 
verkehrenden Elektronenwellen greifen in den Mechanismus der elektrischen Leitfahigkeit, der 
Warmeleitfahigkeit und anderer Transportphanomene entscheidend ein. Da sich jedoch das 
Mikropotential ungeachtet der genannten Wirksamkeit seiner unmittelbaren Vermessung mit 
makroskopischen Geraten entzieht, werden wir uns fortan von seiner Feinstruktur emanzipieren, 
indem wir es durch einen jeweils passend gewahlten, raumlichen Mittelwert ersetzen; dieser ist 
als solcher dem gesamten Zentralgebiet zuzuschreiben. 


b) Unter den Voraussetzungen (11) definieren die Angaben 


low 1 | a 
1<j= <A, LSet —|s- (13) 


J 


oe 
d 


das Randgebiet des metallischen Blockes. Obwohl auch hier, von Sekundareffekten abgesehen, 
das Raumgitter der Atomkerne nicht wesentlich von jenem des Zentralgebietes abweicht, bildet 
sich doch eine deutlich unperiodische Verteilung des elektrischen Skalarpotentiales aus; sie 
wird in einer merklichen Differenz des auf den Grenzflachen des Quaders gemittelten Potentiales 
gegen dessen mittleren Zentralwert manifest, welche ihrerseits, im Gegensatz zu den periodischen 
Mikropotentialschwankungen im Zentralgebiet, der makroskopischen Beobachtung zuganglich 
sind. Um diesem Unterschied auch sprachlichen Ausdruck zu verleihen, mag die Gesamtheit 
dieser beobachtbaren Potentialdifferenzen im funktionalen Begriff des Makropotentiales 
zusammengefaBt werden. 


1* 
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Auf Grund der vorstehenden, vorwiegend deskriptiven Analyse kann keine scharfe Grenze 
zwischen dem Zentralgebiet und dem Randgebiet des Metallblockes angegeben werden. Um 
diesen methodischen Mangel der bisherigen Darstellung grundsatzlich zu beseitigen, vereinigen 
wir die beiden vordem gedanklich getrennten Gebiete zu einem einheitlichen Bereich, ftir welchen 
wir die Giiltigkeit der halbklassischen Elektronentheorie SOMMERFELDs in Anspruch nehmen. 
Zu diesem Zwecke vertauschen wir die Konzeption der Blochschen Elektronenwellen mit deren 
komplementirem Bilde: Der Existenz unzerstérbarer, korpuskularer Elektronen je der invarian- 
ten Ladung (— qo) und der Ruhmasse mp). Wahrend der sozusagen organischen Entstehung des 
kontrollierten Kristalles mége nun immer eines seiner Atome nach dem andern zuwachsen, deren 
jedes, entsprechend seiner OrdnungszahlZ im MENDELJEFFschen System der chemischen Elemente, 
urspriinglich genau Z Elektronen in seiner den Kern umgebenden Hiille enthalte. Vermége der 
angenommenen Genetik lassen sich die Metallelektronen — in deutlichem Gegensatz zu den 
formalen Aussagen der BLOCH-SCHROEDINGERschen Wellenanalyse — gedanklich in zwei Gruppen 
physikalisch wesentlich verschiedener Eigenschaften klassifizieren : 


1. Die dem Kern deseinzelnen Atomes in dessen Urzustand nahmen Elektronen bleiben ihrem 
Mutterkern auch im Kristall dauernd verbunden. 

2. Jedes Atom enthalt in seinem Ausgangszustande 1 < Z’ < Z kernferne ,,Leuchtelektro- 
nen‘, welche als solche durch relativ nur schwache Krafte an den Kern gefesselt sind. Sie ver- 
mégen sich daher nach dem Zusammentritt der Atome zum Metallverband mit Leichtigkeit je 
von ihrem Mutteratom zu lésen. Die Gesamtheit der nach diesem Dissoziationsvorgange ver- 
bleibenden Atomriimpfe bildet ein Raumgitter positiver Ionen der ,,eingeprag.en“ Konzentration 


a 
=a (14) 
deren jedes die invariante Ladung 
Q= HZ (15) 


beherbergt. Dagegen durchschwarmen die freien Elektronen ihren gesamten Lebensraum inner- 
halb und auBerhalb des Metalles als Gas der Konzentration 1, welche sich im allgemeinen von 
Ort zu Ort andert; welche Gesetze regeln ihre Verteilung ? 

Da wir uns in der vorliegenden Arbeit durchaus auf den statischen Zustand des Systemes 
beschranken, mu8 innerhalb des Gases der freien Elektronen Gleichgewicht herrschen. Der 
thermodynamische Aspekt dieser Bedingung wird durch die Annahme einer einheitlichen 
absoluten Temperatur 

i =const. (16) 


der gesamten Gassphare erfiillt. Damit iiberdies das mechanische Gleichgewicht, deran einem 
beliebig gewahlten Gaselement angreifenden Krafte gewahrt werde, mu8 der zufolge (16) iso- 
therme Gradient der GiBBsschen Elektronengas-Funktion y = »(x, y, 2) innerhalb der Gassphare 
verschwinden 


grady =o. (17) 
Diese Gleichung im Verein mit der aus ihr durch Integration hervorgehenden Aussage 
Wes Viz) =p, == Const, (18) 


tritt somit in der an SOMMERFELDs Elektronentheorie der Metalle anschlieBenden, makroskopi- 
schen Elektrostatik dieser Stoffe an Stelle der als unzureichend erkannten Beschreibungen (6) und 
(7) des namlichen Feldes durch die klassische Elektrodynamik. Durch diese tiefgreifende Einsicht 
wird daher das elektrostatische Skalarpotential seiner klassisch beherrschenden Stellung als 
Festwert entkleidet ; diese kommt vielmehr gema8 (18) der Grppsschen Funktion py zu, wahrend 


P = 9%, ¥, 2) (19) 
als vorerst unbekannte Ortsfunktion aufzufassen ist. ; 


Wir richten nun unser Augenmerk auf ein individuelles Elektron, dessen ,, Lebensiauf“ wir die 
NeEwrTonsche Dynamik materieller Punkte zu Grunde legen. Da wir uns in ihr von den HEISEN- 


" 
43 
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BERGschen Ungenauigkeitsrelationen bewuBt emanzipieren, ist es physikalisch sinnvoll, jenem 
Elektron im Kontrollaugenblick ¢ gleichzeitig die drei Ortskoordinaten (x; y; z) und die drei be- 
ziehentlich achsenparallelen Komponenten (v,; v,; v,) seiner resultierenden Korpuskular- 
geschwindigkeit 

v=sor+u+u; (20) 
zuzuschreiben. Die Verteilung dieser Geschwindigkeitskomponenten auf die Gesamtheit der 
freien Elektronen gehorcht, in dem vorausgesetzten Zustande des thermodynamisch-mechanischen 
Gleichgewichtes, der ,,spinmodifizierten“ Fermistatistik: Seien 4 die PLANcKsche und k die 
BoLtzmMaNNsche Konstante, so halten sich im infinitesimal kleinen Elemente dv, dv, dv, des 
sogenannten Geschwindigkeitsraumes 


3 ) 
dn =2 ("") . ee (21) 
; pee oe 
Ex f 
p or + 


freie Elektronen je Einheit des Konfigurationsraumes auf. Zufolge unserer Zustimmung zur un- 
eingeschrankten Benutzung der NEwTonschen Mechanik sind nun bei invariant zu denkender 
Masse m = mp jedes Einzelelektrons fiir jede der drei unterschiedlichen Geschwindigkeitskompo- 
nenten alle reellen Werte zulassig. Daher resultiert aus (21) im Verein mit (16), (18) und (20) fiir 
die Konzentration 

n= n(x, y,2) (22) 
die Integraldarstellung 


Re a ("2 dug dvy dv, 
h 1 2 2 2 
z Mo (vz + vy + V2) — (Po + % Y) 
Exp - cine 
—-wo —wo —oo 


Re 
m, \3 _- v? du 
=2("") 42 ea eS (23) 
| 7. Mg V0 — (Wo + Jo Y) 
J Exp eo oT: +4 


Als quasineutral bezeichnen wir jene Gebiete des Metalles, in welchen bei der Mittelung iiber 
ein zwar ,,physikalisch“ hinreichend groBes, gleichzeitig ,,mathematisch“ aber geniigend kleines 
Raumelement der Erwartungswert der algebraisch aufsummierten elektrischen Ladung verschwin- 
det. Die dort herrschende, , ,eingepragte‘‘ Konzentration 1, der freien Elektronen berechnet sich 
somit aus der Bilanz 

— qn, +QN =0 (24) 
unter Berufung auf (14) und (15) zu 
{LGD 


We Le a 


(25) 

Bei allen technisch realisierbaren Temperaturen T fallt nun diese Konzentration so hoch aus, 
daB sich das sie bildende Gas des freien Metallelektronen im Zustande iiberaus starker Entartung 
befindet. Auf Grund dieses Sachverhaltes gelangen wir zu einer dort gewiB ausreichend genauen 
Beschreibung des gesuchten elektrostatischen Feldes, indem wir durch den allerdings nur ideell 


méglichen ProzeB 
T->o (26) 


zum vollstandig entarteten Elektronengas iibergehen. Obwohl dieses Gedankenexperiment, 
seiner Einfiihrung gema8, zunadchst auf die quasineutralen Gebiete des Metalles beschrankt ist, 
werden wir sie fortan auf den gesamten Block 


|% | 


Deo << [ols 
rah ees Peers 7 (27) 
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ausdehnen; dagegen behalten wir uns vor, ungeachtet der tatsachlich stets wirksamen Gleich- 
gewichtsbedingung (16), in der Umgebung des jeweils kontrollierten Metalles entweder die 
operative Aussage (26) zu benutzen oder aber an der Voraussetzung eines endlichen Wertes 
T > der dort herrschenden absoluten Temperatur festzuhalten. 

Zur Fermistatistik (21) zuriickkehrend, erschlieBen wir aus (26) im Verein mit (20), daB keines 
der freien Elektronen innerhalb des Metalles die Grenzgeschwindigkeit 


Vee = | Pot oP (28) 
Tq 
iiberschreiten kann. Dort reduziert sich daher (23) auf die Angabe 
sitet g ares & ea (tel sie Soe 
n= a() 4% [ av = an 2 a : (29) 


Da nun entweder in der GrBBsschen Funkton oder im elektrischen Skalarpotential m eine 
additive Konstante physikalisch belanglos bleibt, diirfen wirin der Summe (y + qo ~) = (Wo + %) 
fortan ohne Beschrankung der Allgemeinheit 


Vonea0 : (30) 
setzen. Zufolge dieser Wahl vereinfacht sich (29) zu der Relation 
= 8a fm? Ge \*" 
ae (") (2 m9 *) (31) 
mit deren Hilfe wir der eingepragten Konzentration n, nach (25) das eingepragte Potential 9, 
zuordnen: 
__ 8% /m\?(, d a EE ans ek ee 5 3 eo 
Up 2 ("| (29, oe 43 8 h2 | , Pe 8 mM I i é n,) : (32) 


Fiihren wir also n, und ¢, beziehentlich als ,,natiirliche Einheiten“ der Konzentration und des 
Potentiales g im Metalle ein, so lassen sich die Gl.(31) und (32) zu den dimensionsfreien, fiir alle 
Metalle einheitlich geltenden Aussagen 

n — 3/2 : Pp n \2/3 

ae ea (33) 
zusammentassen, welche durch Bild 2 veranschaulicht werden. 

Wir vergleichen sie mit den entsprechenden Eigenschaften eines ,,idealen“ Gases freier Elek- 

tronen, deren Geschwindigkeitsverteilung also durch die Formel 


1 
ae My (Ve + Vy + v2) — (Wo + 9) 
— —| dv, do, de, 


__ , (1M \* 
dy == 2 (“*) Exp AT 


-. My V® — (Wy + Jo) 
AT 42 v2 dv (34) 


geschildert wird; aus ihr folgt durch Integration 
Yo + We 
kT 


my \3 {22k T \8? 
n= 2 °é 
sl Mg 


Stimmen wir abermals der Wahl (30) fiir den konstanten Wert der GrpBsschen Funktion zu, 
messen jedoch diesmal die Konzentration in der ,,natiirlichen Einheit‘ 


(35) 


my \8 [2k T \8? 
ie 2 ) ( ws - el (36) 
und des Potential als vielfaches der ,,Temperaturspannung“ 
kT 
s Ric ss (37) 


i hw Se eee FeE—E Ee aii 
= ———— — oe ~ > — 


~ 
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so entsteht aus (35) die dimensionslose Relation 


NT 


welche in Bild 3 graphisch dargestellt wurde. 


3 
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Bild 2. Zusammenhang zwischen Elektronenkonzentration Bild 3. Vergleich der Konzentrationen im entarteten 

und Potential. und im idealen Gas freier Elektronen. 


3. Beweglichkeit und Diffusion der freien Elektronen 


Im vorstehenden haben wir den Zusammenhang zwischen der Konzentration 2 des Elek- 
tronengases und dem jeweils am gleichen Orte herrschenden elektrischen Skalarpotentiale 
durch Integration der von Fall zu Fall zustandigen Statistik auf wesentlich nur formalem 
Wege hergestellt; der stationare Zustand des Systems wurde hierbei durch die Forderung 
seiner einheitlichen Temperatur im Verein mit dem Gleichgewicht der mechanischen Krdafte 
[verschwindender Gradient der GipBsschen Funktion} verbiirgt. Dieser wesentlich statischen 
Auffassung laBt sich eine dynamische gegeniiberstellen: Wir deuten den stationaéren Zustand 
des Systems als Resultat zweier verborgener kinetischer Vorgange, welche einander gegen- 
seitig kompensieren: Die von den Kraften des elektrischen Feldes erzwungene Konvektions- 
stromdichte 7 werde von der Diffusionsstromdichte s aufgehoben, welche ihrerseits von 6rt- 
lichen Konzentrationsunterschieden herriihre. 

Das in diesem Satze enthaltene Programm ist von EINSTEIN fiir den Fall idealer Gase 
durchgefiihrt worden: 

Es sei f die [positive] Beweglichkeit der Elektronen, D ihr Diffusionskoeffizient im Exi- 
stenzgebiete ihres idealen Gaszustandes. Mit Hilfe dieser Zahlen ergibt sich fiir die Konvek- 
tionsstromdichte 7 die Angabe 


G=QnPBE=—Qqnf grado, (39) 
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wahrend die Diffusionsstromdichte s aus 


S = q D grad n (40) 
zu berechnen ist. Demnach lautet die angekiindigte, dynamische Stationaritatsbedingung 
j+s=—qnfgrady +qDgradn =o, (41) 

oder : 
rad Inn = * grad 2 
grad In m = — grad p . (42) 


Andererseits entnehmen wir aus (35) die Relation 


grad Inn = mc grad, (43) 


deren Vergleich mit (42) das wohlbekannte EINSTEINsche Gesetz 


== (44) 


liefert. 

Um diese Uberlegungen sinngemé8 auf das vollstandig entartete Elektronengas iiber- 
tragen zu kénnen, miissen wir auf die gaskinetische Wurzel der Begriffe Beweglichkeit und 
Diffusion zuriickgehen. 

Wir gehen von dem Fall verschwindenden Makrofeldes aus und weisen dann den freien 
Elektronen unseres Metallmodelles je eine ,,freie Weglange“ / zu, nach deren Durchfliegen 
sie immer mit einem der positiven Ionen elastisch zusammenstoBen médgen; von der mit 
diesem StoBe genetisch verbundenen elektromagnetischen Strahlung sehen wir geflissentlich 
ab. Durch </> die mittlere freie Weglange bezeichnend, nehmen wir an, daB die Wahr- 
scheinlichkeit dw(/) einer im infinitesimal schmalen Intervall zwischen / und (/ + dl) ge- 
legenen freien Weglange dem CLausiusschen Verteilungsgesetz 


aul) =o oH (45) 


gehorche; aus ihm folgt fiir das mittlere Quadrat der freien Weglange die Relation 
= (Pdw 2 Oe, (46) 


von welcher wir spater Gebrauch machen werden. 

Wir denken uns jetzt dem kontrollierten System ein zeitfreies elektrisches Feld der vek- 
toriellen Starke E aufgezwungen, deren absoluter Betrag indessen als so klein vorausgesetzt 
werde, daB sie die schon vordem bestehende Schwarmbewegung der freien Elektronen nicht 
merklich zu st6ren vermag. Eines dieser Elektronen, welches im Beobachtungszeitpunkt ¢ 
die freie Weglange mit der Schwarmgeschwindigkeit v durcheilt, unterliegt somit wahrend der 
Zeitspanne 

ieee (47) 
der Beschleunigung 
Jo 
as es Jos (48) 
so da sich nunmehr der vordem vollig wirren Schwarmbewegung der Gesamtheit aller freien 
Elektronen ein antiparallel zur Feldrichtung ziehender Trift ttberlagert. Nehmen wir jetzt 
an, daB die ihn beschreibende, gerichtete Zusatzgeschwindigkeit v beim ZusammenstoB des 
jeweils kontrollierten Elektrons mit einem der positiven Ionen zur Ganze vernichtet wird, 
so mibt 


ty = asta eS, (49) 
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jene Zusatzgeschwindigkeit am Ende des beobachteten freien Fluges, wahrend 


Oe ete 
2M v 


[~) bo 


an = — , (59) 
die gleichzeitige Verriickung des Elektrons antiparallel zum Feldvektor beschreibt. Aus (47) 
und (50) resultiert fiir den Erwartungswert <v> der Triftgeschwindigkeit der Aus- 
druck 


ee (52) 
/ 


Gah eDrY/ (52) 


vereinfacht. Wie aus seiner Herleitung ersichtlich ist, bezeichnen (3) und (=\ be- 
ziehentlich den Erwartungswert der reziproken quadratischen Schwarmgeschwindigkeit und 
der reziproken Schwarmgeschwindigkeit bei deren Mittelung iiber alle zeitlich nachein- 
ander erfolgenden St6Be eines einzelnen, individuellen Elektrons; als solcher ist er von 
dem Scharmittelwert der reziproken Schwarmgeschwindigkeit bei deren Mittelung iiber 
alle gleichzeitig stattfindenden StoBe der verschiedenen Elektronen gedanklich scharf 
zu trennen. 

Um die beiden, in (52) eingehenden Erwartungswerte explizit zu berechnen, registrieren 
wir die Anzahl M der Zusammenst6Be, die ein individuelles Elektron der Schwarmgeschwindig- 
keit v wahrend der Zeitspanne T —> oo erleidet; sie ergibt sich aus der Definition der mittleren 


- freien Weglange </> mittels der kinematischen Bilanz 


M-(ly=T-v. (53) 


Unter den jeweils freien Elektronen, welche sich je Raumeinheit in der Umgebung des 
Aufpunktes (x, y, z) aufhalten, mége nun der Bruchteil 


f(v) dv, (54) 


eine Schwarmgeschwindigkeit des zwischen v und (v + dv) gelegenen, infinitesimal schmalen 
Intervalles dv besitzen. Da dann die Gruppe dm dieser ,,v-Elektronen‘‘ wahrend der Zeit- 
spanne T entsprechend (53) und (54) gerade 


dn 
nn 


Mn. flo) dv =~ -n-9- flv) dv, (55) 
ZusammenstéBe erlebt, miBt das Integral 
[ Mom: fe) dv= fom [ vf as, (56) 
O° ; O° 


deren Gesamtzahl. Demnach resultieren fiir die gesuchten Erwartungswerte beziehentlich 
die Berechnungsvorschriften 


c (57) 
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und 


so daB (52) in die Formel 


| = flv) do 
VU 


=— BED “= 
& S flv) dv 


(59) 


iibergeht; wir wenden sie auf zwei Sonderfalle an: 


1. Falls sich das Kollektiv der freien Elektronen wie ein ideales Gas verhilt, gilt nach 
(34), (35) und (54) 


= Mo 3/2 ° 2 v= > me 60 
fo) = ("2 42 v € ; (60) 
Daher liefert (59) die Angabe 
Se Se ae 61 
Chinas BD ale (61) 
welcher man die Beweglichkeit 
See a8 ae 6 
Ce ar Ger 2nkT (62) 


entnimmt. 


2. Befindet sich das freie Elektronengas im Zustande der vollstandigen Entartung, 
so folgt aus (21), (26), (28), (29), (30) und (54) die Verteilungsfunktion 


ONS =; = 30 (me) (63) 


er 24P 
Ihre Substitution in (59) fiihrt im Verein mit (31) auf die Relation 


=> 1 
By =—- BEA —2__ = #3 (Z)PE, (64) 
Mo 2 % h 3n 
\: Mo ? 
Die Beweglichkeit \ 
za <u» eal: a \4/3 
i Eh a0) (, 4 ; (65) 


erweist sich nunmehr als Funktion der Konzentration 1, so daB sie gleich dieser im 
allgemeinen von Ort zu Ort verschieden ausfallt; im Gegensatz zur ,,klassischen‘‘ Beweglich- 
keit (62) hangt dagegen die Beweglichkeit im entarteten Gas nicht von der Ruhmasse 
der Elektronen ab. 

Die einfachen Uberlegungen, welche uns zur Kenntnis der Beweglichkeit f$ verhalfen, 
versagen bei der Berechnung des Diffusionskoeffizienten D: Der im makrofeldfreien Zustand 
des Systemes ja typisch irreversible Vorgang der Diffusion ist mit der Voraussetzung 
des thermodynamischen Gleichgewichtes schlechterdings unvereinbar. Erst die kinetische 
Auffassung des Gleichgewichtszustandes erméglicht es, die angezeigte Schwierigkeit auf 
gedanklich einwandfreiem Wege zu iiberwinden: Zur EINsTEINschen Relation (44) zuriick- 
kehrend, finden wir mit Riicksicht auf (62) fiir den Diffusionskoeffizienten D des idealen 
Gaszustandes der freien Elektronen zunachst die Formel 


Da" pah-aVre, (66) 
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welche einer einfachen kinematischen Deutung fahig ist: Nach (54) und (60) miBt 


s aoe ge 


<v» =| v f(v) dv = ES ? (67) 


den Scharmittelwert der in einem festen Kontrollzeitpunkt des Kollektiv der freien Elektronen 
kennzeichnenden Schwarmgeschwindigkeit, so da} (66) in 


D = —- d+ Wy (68) 


iibergeht. 

Wir behaupten nun, da8 die naémliche Relation auch den Diffusionskoeffizienten D des 
Kollektivs freier Elektronen in dessen vollstandig entartetem Zustande angibt, sofern man 
nur die in ihm herrschende Diffusionsstromdichte s, in sinngemaBer Verallgemeinerung der 
klassischen Gleichung (40), durch 

$ = q grad (Dn) (69) 
beschreibt. 

Zum Beweise dieses Satzes bilden wir zunachst aus (63) den Scharmittelwert <v> der 
Schwarmgeschwindigkeit 


Ver 


i» = [ vf) = 2 \ ako (70) 


ie) 


und erhalten gem&B (68), bei erneuter Verwendung des Zusammenhanges (31), die Angabe 


D= =. <ld + <vy ==. <1) \ me? = = Ox ~(2n)". (71) 


m 7 


Wie vordem die Beweglichkeit 6 nach (65), so erweist sich jetzt auch der Diffusions- ~ 
koeffizient D gleichzeitig mit der Konzentration » im allgemeinen als Funktion des 
Ortes. Aus dem Vergleich dieser beiden gaskinetischen Kenngré8en resultiert die einfache 


Relation 
Ds ee h? 3 le eet 
eet a (72) 


welche, wenn man so will, als Ersatz des auf ideale Gase beschrankten EINSTEINschen Gesetzes 
im Falle ihrer vollstandigen Entartung angesehen werden kann. 

Auf Grund der in (59) gegebenen Darstellung der Diffusionsstromdichte s lautet die kineti- 
sche Gleichgewichtsbedingung des vollstandig entarteten Elektronengases 


jts=—nqQf grad y + gy grad (Dn) = 0, (73) 
welche mit Riicksicht auf (65) und (71) in die partielle Differentialgleichung 
3 h BN 
seine) 


ubergeht. Zu ihrer Lésung machen wir nach Wahl zweier vorerst unbekannter Konstanten 
A und a den Ansatz 


1/3 
n > Bees: cS grad m = grad (74) 


g=A (2 nm) (75) 


nh 


und erhalten aus (74) die Forderung 
4, ae een 8) ea ret ie. 4 (ig oad A 
ae dy A «(3 m) grad (3 n) = “ly air (2 n) grad (3 n) ; (76) 


welche mittels der Angaben 


(77) 
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identisch in » befriedigt wird. Demnach stiftet Gl. (75) zwischen dem elektrischen Skalar- 
potential m und der Konzentration des vollstaindig entarteten Elektronengases in dessen 
Gleichgewichtszustande den Zusammenhang 


o=z_ (2a), (78) 


welcher inhaltlich genau mit der Aussage (31) der FERMI-Statistik iibereinstimmt; gerade 
diese Gleichheit enthalt den gewiinschten Beweis der Gl. (71). Indessen gestattet es der 
vorstehend entwickelte Gedankengang, iiber die notwendig nur zeitfreien Aussagen des sta- 
tistischen Gleichgewichtes prinzipiell hinauszugehen: Fassen wir 6 nach (65) und D nach (71) 
je als Funktion der Konzentration m auf, so vermag die partielle Differentialgleichung der 
Diffusion 

én 


<p = V2 (D+ n) — div (n B grad g) , (79) 


im Verein mit der von der elektrischen Raumladungsdichte @ diktierten PorssoNschen 
Potentialgleichung auch dynamische Vorgange beliebigen zeitlichen Verlaufes im voll- 
standig entarteten Gase der freien Elektronen zu beschreiben; allerdings darf nicht ver- 
schwiegen werden, da wir uns bei einer solchen Erweiterung einer gewissen Inkonsequenz 
schuldig machen: Der Begriff der Grenzgeschwindigkeit v,, nach (28), welche wesentlich 
in die Berechnung von f und D eingreift, entstammt ja der im Grenzfalle verschwindend 
niedriger absoluter Temperatur resultierenden Fermistatistik, die ihrerseits nur fiir den 
Zustand des Gleichgewichtes Giiltigkeit beansprucht! Da wir uns jedoch im Rahmen der 
vorliegenden Arbeit durchaus auf stationaére Vorgange beschranken, ist die angezeigte Dia- 
lektik fiir unsere Zwecke irrelevant. 


4. Das elektrostatische Feld im eben begrenzten Halbkérper 


Durch den GrenzprozeB 
A>OO; UO; YOu, (80) 


gehen wir zu einem ideellen, allseitig unbegrenzten Metallblock iiber, welcher langs der Ebene 

z =o durchschnitten werde. Denken wir uns jetzt den Halbk6érper z < 0 beseitigt, so ent- 

wickelt sich in dem verbleibenden Kristall z > 0 ein elektrostatisches Feld, dessen Potential » 

wir aufzusuchen haben. Bezeichnen wir durch A, die sogenannte absolute Dielektrizitats- 

konstante des leeren Raumes und durch ¢ die relative Dielektrizitatskonstante des Metalles, 

so lautet die Porssonsche Gleichung 

Pie es 

Vp ee (81) 

In ihr resultiert @ aus der ,,makroskopischen“ Dichte @, der positiven Ladung, welche 
entsprechend (24) ein fiir allemal durch das Produkt 


0, =ON=Qn, (82) 
gegeben ist, und der Dichte o_ der negativen Ladung, welche sich ihrerseits nach MaBgabe 
der jeweiligen Konzentration » der freien Elektronen im allgemeinen von Ort zu Ort andert 


0- = —Jn- (83) 
Stellen wir nun der in (32) definierten Potentialeinheit m, durch 


Age P 
NOS 6 
an ea (84) 
das Quadrat der natiirlichen Langeneinheit L im Kristall zur Seite, so geht (81) in die dimen- 
sionsfreie Gestalt 
My es ela ea 
eA eel 5) 


iiber. 
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Der Einfachheit halber beschranken wir uns weiterhin auf ebene Potentialfelder, welche 
als solche nur von der z-Koordinate abhangen. Nach Einfiihrung der ,numerischen‘ 
Koordinate 


iar (86) 
und des ,numerischen Potentiales“ 
AEE at (87) 
reduziert sich dann (85) auf die gewohnliche Differentialgleichung 
coed 1 | 
eee ae ed 


bei deren Lésung wir drei Falle zu unterscheiden haben: 


1. Im Kontrollgebiet herrscht ElektronentiberschuB 


Sah (89) 


nN 6 


2. Im Kontrollgebiet herrscht Elektronenmangel 


W 
Ones ae (90) 


é 


3. Die freien Elektronen sind aus dem Kontrollgebiet véllig verdraéngt worden 
ee (91) 


Ungeachtet der unterschiedlichen Voraussetzungen (89) und (go) resultiert in beiden 
Fallen aus (87) im Verein mit (33) fir @ die gleiche, nichtlineare Differentialgleichung 


a® 1 . 
pea a tly (92) 


wahrend im Falle (go) das numerische Potential ® der linearen Differentialgleichung 


ad 1 
Soe eee (93) 
unterliegt. 
Die Behandlung der Gl. (93) zuriickstellend, erganzen wir (92) durch die Randbedingungen 
fim <= 0 (94) 
o> O% 
und 
lim 6 = @ = (“° 6 
to Tig ee ea : z (95) 
in deren letztgenannter also 
ne en), (96) 
¢—->o 


den Grenzwert der freien Elektronenkonzentration m an der Front € = o des Kristalles angibt, 
falls man sich dieser von innen her nahert. 
Unter Berufung auf die Identitat 
oS = ne E @ 
ac? d®@ 


d® 


ae 


; B= 


Pe 


liefert (92) mit Riicksicht auf (94) durch einmalige Integration 
ee Pf P). (98) 
ro) 


— 
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Von hier aus gabelt sich der Gedankengang: 
1. In Gebiete des Elektroneniiberschusses - > 1 gilt gem4B (33), (87) und (89) auch ® > 1, 
WL 


sodaB dort nach Ausweis des Bildes 4 die in (98) definierte Funktion /(®) stets positiv aus- 
fallt. Soll dann das elektrostatische Feld im Halbkérper ¢ > 0 beschraénkt bleiben — und 
dies ist von einer physikalisch sinnvollen Lésung gewiB zu verlangen — so muB also das 


2 


45 


0 0,5 4 4,5 g an) 


o— 


Bild 4. Das Quadrat des numerischen Feldgradienten als Funktion des numerischen Potentiales. 


numerische Potential ® von seinem Randwerte ®) > 1 an nach dem Innern des Metalles hin 
monoton absinken. Im Einklang mit diesem Schlusse finden wir aus (98) nach Trennung der 
Veranderlichen das Integral 


(= | —— =. (99) 


welches allerdings nur mittels numerischer Berechnung explizit hergestellt werden kann. 
Das Ergebnis wird durch Bild 5 veranschaulicht, in welchem der Gang ® = @(¢) des numeri- 
schen Potentiales fiir unterschiedliche Randwerte ®) graphisch dargestellt ist; in jedem Falle 
nahert es sich ,,von oben her“‘ dem Grenzwerte 

lim Dies m; (100) 


f>o 


so daB ebendort der Kristall elektrisch quasineutral ist. 
2. Im Gebiete o< <1 des Elektronenmangels unterliegt auch das numerische Po- 
nN 


tential ® der Ungleichung o << ® < 1, so dal} wir an Hand des Bildes abermals auf den gleich- 
zeitig positiv-definiten Wert der Funktion /(®) nach (98) schlieBen. Der zusatzlichen Forde- 
rung einer im Halbkérper ¢ > o beschrankt bleibenden Lésung mit EinschluB ihrer ebendort 
zu wahrenden Stetigkeit geniigen wir dann nur mittels einer Funktion ®, welche mit zu- 
nehmendem € monoton ansteigt. Demnach folgt aus (98) das Integral 


Pp 
dD 
seg Se ea (101) 
3 235 D 
ps r7D) 
ee 
2, 
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Seine numerische Berechnung liefert bei unterschiedlich vorgegebenen Randwerten 
o < ®) <1 den raumlichen Gang ® = @(¢) des numerischen Potentiales, welches sich ent- 
sprechend Bild 5 ,,von unten her“ abermals dem Grenzwerte 

lim @=1 (102) 
¢—>0o 
nahert; wiederum erweist sich dort das Metall als quasineutral. 

Nach den Ergebnissen der vorstehenden Uberlegungen fallt sowohl bei UberschuB freier 
Elektronen im Metalle wie bei deren Mangel der Randwert ®, des numerischen Potentiales 
stets positiv aus. Um daher im Metalle elektro- 
statische Felder der Eigenschaft ®, < o reali- 
sieren zu kénnen, haben wir eine Randschicht 
der zunachst noch unbekannten [numerischen] 
Dicke 6 zuzulassen, welche von freien Elektronen 
vollig entbl6Bt ist. Fiir das numerische Poten- 
tial des Bereiches 0 <p <0 ist somit die Diffe- 
rentialgleichung (93) zustandig, wahrend das 
Feld in¢ = 6 nach wie vor der Differentialglei- 
chung (92) unterliegt. Zufolge (33) konvergiert 
nun das numerische Potential ® mit gegen 6 
abnehmendem ¢ gegen ® = 0; daher findet man 
das in € = 6 herrschende Potential ® aus (101), 
indem man ® =o setzt und iiberdies den 
Ursprung des Bezugssystemes um 6 verschiebt: 


DP 
7D 
a0 Mie -. (103) 
SS) 2? Gil2_@ 
le 


d@ 
dc 


Fin 6 <¢€— 0 konvergiert somit der Gradient 
gegen den Grenzwert 


—s 0 2 5 78 410 
> olen ==> 
(Pps cali (104) is 
C6 G 5 Bild 5. Verlauf des numerischen Makropotentials im Metall. 


In der Ebene € = 6 mu8 nun sowohl das in € < 6 herrschende numerische Potential ® 
wie dessen Gradient a. stetig in die entsprechenden GréBen des Bereiches ¢ = 6 iibergehen. 


Wir geniigen dieser Forderung durch das auf die Randschicht beschrankte Integral 


es ere reaped le (6—¢), (105) 
4 a 
der Differentialgleichung (93), welchem wir den Randwert 
%=— 8-3 d<0, (106) 
A 5 


des numerischen Potentiales entnehmen; ist dieses vorgegeben, so ergibt sich aus (106) um- 
gekehrt die Dicke 6 der Randschicht zu 


Sree 5 
i= 2/3 / so i> ) (107) 


entsprechend Bild 6. Aus (105) resultiert der Gradient 
d@ ee! 30 
dt mes (0 ¢) I Ie ? (108) 
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welcher in der Ebene € = 0 den Wert 


dD 1 B 3 5 
fel-ze+ | -\/3 \): ye (209) 


annimmt; man entnimmt dieser Gleichung die einfache Relation 


ee Spee (110) 
oO 5 5) 


at j 


= Sa 

-2,9 =A =r) =| =10)5) 0 

eres 

Bild 6. Numerische Dicke der Randschicht als Funktion des (negativen) numerischen 
Randpotentiales im Falle der Elektronenerschopfung. 


5. Numerische Angaben 


Um uns mit dem quantitativen Inhalt der bisher erzielten Ergebnisse vertraut zu machen, 
mégen die in ihnen auftretenden, sozusagen natiirlichen Einheiten fiir das Beispiel chemisch 
reinen Kupfers zahlenmaBig berechnet werden; sein Atomgewicht wird, in hier hinreichender 
Genauigkeit, durch 

A = 63,8 
gemessen, wahrend 
kg 
Vy ==8,8 « 10° —~ 
seine Dichte angebe. 


Aus (9) folgt nun zunachst im Verein mit (8) der Abstand dzweier Nachbaratome im kubischen 
Modellgitter zu 


3 
638 
i= = = 2,28- 107m 
8,8 - 108 + 0,602 + 107? 4 


so daB gemaB (14) die Zahl 


1 0,602 + 102? + 8,8 + 103 1 
= N=-— = 000s Oe 
a 63,8 m3 


der eingepragten Konzentration der positiv geladenen Atomriimpfe gleicht. Der Einfachheit 
halber werde weiterhin 
Zt 


angenommen, so da gemaf (25) auch die eingepragte Konzentration der freien Elektronen 
durch 


1 
N, = 0,083 - 10% — 
m?® 


eR I pan a a ad HE 
, - 
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' gemessen wird. Mit Hilfe der Werte der PLanckschen Konstanten 


h = 6,63 - 10-* Joule sec 
des absoluten Betrages 
Jo = 1,60 - 107! Coulomb 


der Ladung des Elektrons und seiner Ruhmasse 
Me — 0,4) 105 ke 


berechnet sich dann aus (32) die Potentialeinheit m, zu 


6,63 > 10-88 3 
”= 


8-9,1° 107%! - 1,60 - 10-18 


2/ 
0,0083 - 10%) ps 6,18 V. 


7 


Erganzen wir die Angabe der universellen Konstanten 


A oo =: ee 
durch die allerdings nicht willkiirfreie Annahme der relativen Dielektrizitatskonstanten 
e=2 
so folgt fiir das Quadrat der Langeneinheit L aus (84) 


1+ 6,18 
470°9:* 109+ 2° 1,60° 10719 - 0,083 - 10°09 


= 


= 0,206 - 107-29 m? 


also 
LE =='0;455 7 107 m= 0,2 @5 


6. Uber die Natur der Austrittsarbeit 


Jede in sich geschlossene Theorie der Metalle muB imstande sein, tiber den Mechanismus der 
Elektronenemission Rechenschaft abzulegen. Man geniigt dieser Forderung durch die Annahme, 
daB jedes freie Elektron beim Verlassen seines Muttermetalles eine dieses kennzeichnende Aus- 
trittsarbeit W, zu iiberwinden hat; in der Tat fiihrt diese Hypothese zu einer phanomeno- 
logisch recht befriedigenden Beschreibung sowohl der thermischen wie der lichtelektrischen Elek- 
tronenemission. Um die Genetik dieser Austrittsarbeit zu verstehen, pflegte man das eben 
emittierte Elektron durch ein virtuelles ,,Positron“ zu erganzen, dessen Lage aus jener des 
Elektrons durch Spiegelung an der Grenzflache Vakuum—Metall hervorgeht: die zwischen diesen 
antipolaren Punktladungen wirksame ,,Bildkraft“ ist es, deren langs des Elektronenweges 
erstrecktes Integral die Austrittsarbeit bestimmt. Allerdings verlangt die Festlegung der jeweils 
zustandigen Integrationsgrenzen Zusatzannahmen, welchen nicht die von einer strengen physi- 
kalischen Theorie zu fordernde Uberzeugungskraft innewohnt; ihre Aussagen liefern denn auch 
nur einen gréBenordnungsmaBig brauchbaren Anhalt fiir die Austrittsarbeit. 

Ungeachtet der angezeigten Mangel der quantitativen Durchrechnung kann man doch, die 
Grundannahme der Theorie als richtig supponierend, die Austrittsarbeit W, durch den Vergleich 
gemessener Emissionskennlinien mit daran theoretisch zu erwartender Gestalt experimentell er- 
mitteln. Das Ergebnis pflegt man mittels der Relation 


Wi= 4, (113) 


auf die Austrittsspannung U, umzurechnen, welcher — im Rahmen der klassischen Elektro- 
statik! — eine tiberaus anschauliche Deutung gegeben werden kann.: Entsprechend Bild 7 wird 
die Grenze Metall-Vakuum von einer homogenen elektrischen Doppelschicht eingenommen, 
in welcher das Potential y von seinem Metallwert sprunghaft um das MaB 


Ap = — U, (112) 
auf seinen Vakuumwert abfallt. 


LS) 
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Wahrend nun die Konzeption einer solchen unstetigen Potentialanderung in der klassischen 
Elektrostatik als durchaus legal angesehen werden kann, widerspricht sie den Grundsatzen der 
hier behandelten, elektronentheoretischen Elektrostatik der Metalle; vielmehr tritt an Stelle der 
Doppelschicht und ihres Potentialsprunges der stetige Abfall des Potentiales y von seinem einge- 
pragten Werte gy, in den quasineutralen Gebieten des Metalles auf den Grenzwert y > + 0 im 
Vakuum, welcher zufolge der Relation (33) zwischen dem Potential des Kontrollortes und der 
ebendort’ herrschenden, notwendig positiven 
Konzentration der freien Elektronen nirgends 
“ unterschritten werden kann. Im Lichte dieses 
Metall 72 Vakuum Sachverhaltes sollte man vermeinen, daB die 
Austrittsspannung U, mit dem eingepragten 
Potential w, identisch sei. Dieser zunadchst so 
einleuchtende Schlu8 wird jedoch von der Er- 
fabrung widerlegt: GemaB 

ey (133) 

Bild 7. Schematische Darstellung der hypothetischen Doppel- liegen die gemessenen Werte der Austrittsspan- 
Scene at BN Roath eg Sirah are as 7 Yaver- nung U_, stets merklich hoher als das eingepragte 
Potential y,! Wie kann man diese Ungleichung 

dem physikalischen Verstaéndnis erschlieBen ? 

Seiner Definition nach miBt das eingepragte Potential y, den zeitfreien, réumlichen Mittel- 


wert des in den quasineutralen Gebieten des Kristalles dreifach-periodischen Mikropotentiales 
gy: Durch 


Emissionsrichtung 


dV = dx dy dz (114) 
die jeweiligen Elemente eines dort gewadhlten, hinreichend groBen Kontrollvolumens V bezeich- 
nend, gilt 

w= Z| [ [ora dxayae. (115) 
(V) 


Dieser starren, topographischen Durchschnittshéhe des Potentialgebirges m = (x, y, 2) 
stellen wir nun den seinem inneren Wesen nach mit der Frequenz der Schroedingerwellen rhyth- 
misch durchpulsten, wenngleich auBerlich-formal ebenfalls zeitfreien Erwartungswert <@> 
des Skalarpotentiales y zur Seite, welcher von den ruhelos im Labyrinth des Kristallgitters um- 
herirrenden, freien Metallelektronen ,,erlebt‘‘ wird: Aus der Dichte w = w(x, y z) der Aufent- 
haltswahrscheinlichkeit eines solchen Elektrons im Element dV des oben genannten Kontroll- 
volumens V berechnet sich <g> mittels der Vorschrift 


<o> = SSS (x, y, 2) w(x, y, 2) dx dy dz (116) 
(V) 
im Verein mit der Vollstandigkeits-Relation 
ff fee y, 2) dady dz =. (116a) 
(V) 


Aus dem Aufbau der Atome geht nun hervor, dai gm in der Umgebung der Atomkerne am 
groBten ausfallt. Da iiberdies gerade dort auch die Wahrscheinlichkeitsdichte w ihren Héchst- 
wert erreicht, resultiert der Erwartungswert <p> gréBer als das eingepragte Potential m,. Auf 
Grund dieser Einsicht liegt es nahe, die in (113) definierte Austrittsspannung U, mit <p> zu 
identifizieren 

U,= >. (117) 


Allerdings verandert diese neue, vielleicht gewagte Auffassung die tibliche Konzeption der Aus- 
trittsarbeit in physikalisch recht einschneidender Hinsicht: Fiir die Emissionseines Metall- 
elektrons ist nicht allein die Energetik seines schlieBlichen Austrittes in das Vakuum verant- 
wortlich zu machen, sondern ein wesentlicher Anteil der insgesamt aufzuwendenden ,,Emissions- 
arbeit entfallt auf den ,,inneren‘‘ Ubergang jenes Elektrons von den kernnahen Gebieten des 
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Kristalles in den eigentlichen Lebensraum der freien Elektronen. Ungeachtet aller gewi8 be- 
rechtigten Zweifel an der vorgeschlagenen Deutung (117) des Erwartungswertes <g> kommt ihr 
doch jedenfalls die erzeugende Kraft einer weitreichenden Arbeitshypothese zu; beispielsweise 
gestattet sie, den Besonderheiten der Elektronenemission aus diinnen metallischen Schichten 
quantitativ nachzugehen. — Die explizite Berechnung des Erwartungswertes <g> verlangt die 
Kenntnis aller jeweils fiir die Metallelektronen zustandigen Lésungen ihrer SCHROEDINGER- 
gleichung, welche jedoch nicht in den Rahmen der vorliegenden Untersuchung fallen; demnach 
haben wir uns weiterhin mit der Ungleichung (113) zu begniigen. 


7. Das Kathoden-Grenzfeld 

Die Grenzflache € = o des den Halbraum ¢ > 0 erfiillenden Metalles bilde die negative Elek- 
trode eines Systemes zweier planparalleler Platten, zwischen welchen ein homogenes elektrisches 
Feld der Starke Ey ausgespannt sei. Messen wir diese in der , natiirlichen“ Einheit 
2, 
ee 
welche als ,,eingepragte‘‘ Feldstarke bezeichnet werden mag, so unterliegt die dielektrische 
Induktion an der Grenze des Metalles gegen das Vakuum der Stetigkeitsbedingung 


E,= (118) 


a AD. L |. dp E 
] = — —— ae 
Ave ao a Aes ae 7 Ay =; eS 0. 3 
(119) 
Aus ihr ergibt sich im Verein mit (g8) fiir das 
numerische Randpotential 29 1 
@O(o) = B>1 (120) 


die Gleichung 
3 2 pil2 D dig 2 
le =) : i SOG 0 2 


Nach Ausweis des Bildes 8 nimmt also so- 
wohl das numerische Randpotential ®y selbst | e 


wie die mit ihm durch die Relation = 
lim n = m = 0, BI? > n, (122) 
(oe) 
verbundene Randkonzentration der freien Me- 1 T 


tallelektronen mit steigendem Betrage des Feld- 


starkenverhaltnisses = monoton zu. Diese Orts- 
é 


as : 0,5 + 
beschrankung zieht nun, auf Grund der HEISEN- : 
BERGschen Ungenauigkeitsrelationen, cine mit 
wachsender Randfeldstarke Ey immer krasser 
ausgepragte ,,Platzangst“ jedes Einzelelektrons 
‘ + 4 x = 0 0,5 4 15) 2 
nach sich, welche sich in dessen Maximalge- re 
. . . Cty ae oo 
schwindigkeit Ee 
q z Bild 8. Das numerische Randpotential der Kathode als Funktion 
j oe 0 OQ. der Vakuum-Feldstarke. 
lim. 0) —= |/ 2 — 9 Qo - (123) 
z>0 My 


kinematisch offenbart. Falls daher das Vakuumfeld soweit gesteigert wird, daB das Rand- 
potential y) = y, Dy entsprechend der Ungleichung 

Po = Pe Dy a U, (124) 
die Austrittsspannung tibertrifft, konnen die im Sinne der FERMI-Statistik jeweils energiereichsten 
unter den freien Metallelektronen ihren Wirtsk6rper verlassen: Die Gleichung 


[2 tie 
pe eB | iii [eI ve (125) 


S 5 \Pe Ve 
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definiert gem&B (121) und (124) jene ,,kritische Feldstarke, bei welcher die ,,Feld- 
emission“ der freien Metallelektronen einsetzt. Diese grundsatzlich so wichtige Erscheinung 
bedarf also zu ihrer Erklarung nicht des wellenmechanischen Tunneleffektes durch einen klassisch 
undurchdringlichen Potentialwall nach Bild ga, sondern folgt zwanglos aus der hier durchge- 
fiihrten Integration der Porssonschen Differentialgleichung, welche zu dem durchaus stetigen 
Verlauf des Potentiales nach Bild gb fiihrt! Allerdings darf nicht verschwiegen werden, da8 auch 


Emissionsrichtung 


Metal| —» 


—— Vakuum 


Bild ga. Klassische Vorstellung iiber den Verlauf des Potentiales Bild 9b. Neuer Vorschlag zum Verlauf des Potentiales 
beim Einsetzen der Feldemission. beim Einsetzen der Feldemission (schematisch). 


unsere Theorie sowohl in der Definition (32) des eingepragten Potentiales y, wie in der Definition 

der natiirlichen Langeneinheit L nach (84) vermittels des expliziten Auftretens der PLANCKschen 

Konstanten / einen typisch quantenmechanischen Zug enthalt. Sei beispielsweise Uae 15 
Pe 

so erhalt man mit den frither angegebenen Zahlenwerten ftir Kupfer 


ie Ree IW + — (15) —15]= 6,18 «100 YO 


: 0,455: 10°20 5 m 


Volt 


VAN DER ZIEL? nennt Emissionsfeldstarken von der GréBenordnung 10° , welche, wenn 


gleich merklich niedriger, sich doch nicht allzusehr von unseren Angaben unterscheiden. 

Es bedarf nur einer leichten Erweiterung der vorstehenden Uberlegungen, um fiir Feld- 
starken LE, > E,, den Gang der Emissions-Stromdichte als Funktion eben dieser Feldstarke 
kennen zu lernen; da jedoch der Emissionsvorgang als solcher nicht mehr dem Gebiete der Elektro- 
statik angehort, mag die angedeutete Rechnung einer spateren Arbeit vorbehalten bleiben. 

Zum besseren Verstandnis der hier vorgeschlagenen Auffassung der kritischen Feldstarke 
darf vielleicht daran erinnert werden, daB bei der thermischen Elektronenemission die Polaritat 
der Gliihelektrode nicht von vornherein festliegt, sondern erst von den jeweils vorliegenden Be- 
triebsbedingungen des Systemes diktiert wird: Solange man unterhalb des Sattigungs- 
stromes arbeitet, fungiert die Gliihelektrode ungeachtet ihrer technologischen Bezeichnung 


1 ALDERT VAN DER ZIEL, Solid State Physical Electronics (S. 113) Prentice Hall, Engiewood Cliffs, N.J.: 
1957- 
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als ,,Gliihkathode“ tatsachlich als Anode, deren Randpotential notwendig héher liegt als 
jenes derim Entladungsraum befindlichen , virtuellen‘‘ Kathode [Ebene des Potentialminimums]! 
Erst im Operationsgebiete des ScHoTTKy-Effektes wird die Gliihelektrode zur wahren 
Kathode im physikalischen Sinne dieses Wortes, so daB man von hier aus den gedanklichen 
AnschluB an die Feldemission zu vollziehen hat. 


8. Der ideale Kondensator 


Wir erganzen die im vorigen Abschnitt behandelte Kathode des planparallelen Platten- 
systemes durch dessen Anode, deren Grenzebene gegen das Vakuum den Abstand d gegen jene 
der Kathode aufweise. Solange dann die interelektrodische Feldstarke E der Bedingung 

= eg =e (126) 
geniigt, kann keines der Metallelektronen in das Vakuum iibertreten, so daB in diesem mangels 
elektrischer Ladungen gewiB keine Strémung zustande kommt. Da nunmehr das elektrische 
Feld im gesamten System seinen statischen Charakter wahrt, sind wir berechtigt, es als idealen 
Kondensator zu definieren; bei seiner analytischen Beschreibung werden wir uns der Kiirze 
halber auf ein ,,symmetrisches“ System beschranken, dessen antipolare Elektroden aus dem- 
selben Metall bestehen. 

Unter Berufung auf Gl. (121) diirfen wir das numerische Randpotential der Kathode, welches 
hier mit dem Symbol @, bezeichnet wurde, der Relation 


| a aie (127) 


nach Bild 8 entnehmen. GemaB (126) treffen wir daher an der vakuumseitigen Grenzebene der 
Anode das numerische Potential 


Bee 
Ne a eae (128) 


an. Bei der physikalischen Deutung dieser Formel hat man zu beachten, da zufolge der hier 
gemachten Voraussetzungen Kathode und Anode durch einen ,,nichtanalytischen™ Bereich 
getrennt sind, in welchem die Relation (33) zwischen der Elektronenkonzentration m und dem 
Potential y auBer Kraft gesetzt ist. Daher diirfen wir die fiir das Innere der Anode urspriinglich 


Kathode 


Bild 10. Potentialverlauf im idealen Kondensator (nicht mafistablich). 


benutzte Polentialbasis durch eine willkiirlich wahlbare, neue Basis ersetzen, ohne daB die friiher 
angegebene Feldstruktur eben dieses Gebietes von dieser ja nur formalen Operation beriihrt 
wird. Bezeichnen wir den auf die Einheit g, bezogenen Wert dieser Basis durch das Symbol u 
so miBt also das Produkt 

U=u- Pe (129) 
entsprechend Bild 10 den Grenzwert der Potentialdifferenz zwischen den unmeBbar weit von 
einander entfernten quasineutralen Gebieten beziehentlich der Anode und der Kathode oder, 
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zwar weniger genau, doch kiirzer gesagt, die Kondensatorspannung. Bei der expliziten Be- 
rechnung von u haben wir zwei Falle zu unterscheiden: 


1. Die Vakuumfeldstarke E, sei schwacher als das Kathoden-Grenzfeld E,,, nach (125) und 
geniige iiberdies der Ungleichung 
aoe le (130) 


derzufolge am vakuumseitigen Anodenrand Elektronenmangel herrscht. Unter solchen 
Umstanden garantieren wir die an der Grenzebene Anode-Vakuum zu fordernde Stetigkeit der 
normalen elektrischen Induktion, indem wir die Differenz 


0<@,=0,—u<1 (131) 
unter Berufung auf (98) der Gleichung 
te y a a = ts bei 
1 %o V3 4 2 @, =a —(@,—1) (132) 
bE, 5 5 ; + 


unterwerfen; ihre aus Bild 11 zu entnehmende Lésung liefert im Verein mit (128) die gesuchte | 
numerische Basis w. 


ae Eine | Kathode 
ee 5 


oF 
ee Pale aa 
=; erschépfung 
=o 
Bild 11. Numerische Potentialdifferenz innerhalb der Anode Bild 12 Numerische Zusatzspannung der Elektroden 


als Funktion der Vakuum-Feldstarke, zur Vakuumspannung des idealen Kondensators, 


2. Die Vakuumfeldstarke Ey sei zwar, wie vordem, schwacher als das Kathodengrenzfeld E,,,; 
nichtsdestoweniger mag sie die Ungleichung 


E cs 
E =& Ie (133) 


befriedigen. Die vakuumsceitige Randschicht der Anode ist nunmehr von freien Metallelektronen 
vollig entbl6Bt, so daB wir aus (110) fiir die Differenz 


O, = 9,—u<o (134) 
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die Aussage 


cB 1 Ey 
®,—u=3 (2%) (135) 
oder 
at 1 E,\? 3 
= DY + (2 z)-2 (136) 
gewinnen. 
Um die Ergebnisse (127), (132) und (135) iibersichtlich zusammenzufassen, schreiben wir 
U=9,u=E,-d4+AU (137) 
wobei 
AU = 9, [®, — ®4] (138) 


die zur ,, Vakuumspannung “‘(E: d) additiv hinzutretende, innerelektrodische Spannung angibt ; 
ihre Gr6éBe, gemessen in der Potentialeinheit y,, wird durch Bild 12 als Funktion des Verhialt- 


nisses (-3) graphisch dargestellt. Man erkennt aus dem Diagramm den wesentlich nicht- 
é é 
linearen Charakter der funktionellen Beziehung zwischen der je auf das Innere der beiden anti- 


polaren Elektroden ,,verschmierten‘’ Ladungsdichte 
o=+AE, (139) 

welche man je Einheit der Elektroden-Grenzflachen in Rechnung zu stellen hat (Gau8scher 
Satz!); daher kann man dem hier behandelten, idealen Kondensator keine feste Kapazitat zu- 
schreiben, sondern diese ist von der Gesamtspannung U zwischen den Elektroden abhangig. Un- 
geachtet dieser grundsatzlich wichtigen Einsicht ist allerdings der genannte Effekt bei den iib- 
lichen Vakuumfeldstarken quantitativ so schwach, daB er haufig auBer acht gelassen werden 
darf; erst bei hédchsten Anforcderungen an die Genauigkeit, wie sie beispielsweise an die numerishe 
Berechnung der Eigenfrequenz eines fiir eine Atomuhr bestimmten Hohlresonators gestellt 
werden, mag die angezeigte Abhangigkeit der Kapazitat von der jeweils wirksamen Spannung 
bedeutungsvoll werden. 

Wir erganzen diese allgemeinen Uberlegungen durch die Berechnung der in den Elektroden 
aufgespeicherten Freien Energie. Ausgehend von dem Ansatz 


®=1+/A90 (140) 
des numerischen Kathodenpotentiales [6 > 1] folgt fiir hinreichend kleine Betrage 
AD <1 (141) 
der Differenz A® aus (98) mittels binomischer Entwicklung des Radikanden 
((@)= +20" 6 = AO 4... (142) 
J 


Bei Beschrankung auf ihr vorstehend explizite angegebenes Anfangsglied gehorcht somit ® 
der Differentialgleichung 


cf = —(@—1) V3 (143) 


is 2 
unter der Randbedingung 
D(o) = @,. (144) 
Aus (143) und (144) entnimmt man fiir den Gang des numerischen Kathodenpotentiales die 
approximativ giiltige Darstellung 


les 
Dei Oe) a (145) 
aus welcher man riickwarts den Gradienten 
d® oe oe V3 
Se ae ae (146) 
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berechnet. Je Flacheneinheit beherbergt also die Kathode die Freie Energie 


Lede aes - €Ay ” 9% ic: ase “eds Pe 2 V3 
Me ee (@) dz = —_: 7 (Pu 1) - é dé : 7 (Px 1) ore (147) 


\ 


Gem4B (146) gilt nun am Kathodenrand 


d@ epiat Ey gs Bee 


i = —E D —1)° => — Ol =f 1 8 
Gees ig (9, ) 2 E, Pe (148) 
so daB (147) in die einfache Gestalt 
"Sal (149) 
2€ 13 


gebracht werden kann. 
Um die Freie Energie der Anode zu berechnen, vertauschen wir in (140) das numerische 
Potential © mit der Funktion 


® =G@—@,=14+ 40 <1 (150) 


|AD’| <1 (151) 


vorausgesetzt werde; daher liefert (98) nunmehr durch binomische Entwicklung die Angabe 


in welcher, entsprechend (141), 


f@)=—Aaw 4... (152) 


Sei dann ¢’ der numerische Abstand eines inneren Anodenpunktes von der Grenzebene der Anode 
gegen das Vakuum, so gehorcht ®’ der Differentialgleichung 


do’ ave 
e =o) (153) 
unter der Randbedingung 
P'(0) = D4. (154) 
Ihr Integral lautet somit 
Ae rl V3 
@’ =1—(1— Dy)e a (155) 
so daB die Gleichung 
V3 
dg’ , 3 =H Be 
1 = — 04) 3 (156) 


den raumlichen Gang des Gradienten beschreibt. Demnach enthalt die Anode je Flacheneinheit 
die Freie Energie 


d co 
eA dg \2 eA, pe Bee Weim te chal ane Fee ial ie} 

w= 2 | (yas 2 (aye Sf ee VF ag = 240% Qa. a57) 

Mit Riicksicht auf (156) gilt nun am Anodenrand 

dp’ es Ge Senne 
ipraene CS aise ae ioe) 
so daB (157) in 
A 1D 

Wn a Ey Fy (159) 


umgeschrieben werden kann. Da nun der Interelektrodenraum je Flacheneinheit die Energie 


A 
Wes eo Red (160) 
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Cr 


beherbergt, entfallt auf die Flacheneinheit des gesamten Systemes die Freie Energie 
W = W, + Wo W, =" Bila + 222 

2 € V 3 

In der bisher geiibten Genauigkeit miBt nun gem4B (137), (138), (148) und (158) die Summe 


; (161) 


U=,|a@+—*2| (162) 


die resultierende Kondensatorspannung. Aus (161) und (162) folgt die Relation 


so daB 


2 (164) 
| 


den Grenzwert der Kapazitat je Flacheneinheit des Kondensators im Falle einer infinitesimal 
kleinen Kondensatorspannung definiert; er ist also stets kleiner als der entsprechende Wert, 
welcher aus der klassischen Elektrostatik gefunden wird. 


9. Zusammenfassung 


Die klassische Elektrodynamik beschreibt das elektrische Feld innerhalb eines homogenen und 
isotropen Stoffes mittels nur zweier phanomenologischer Festzahlen: Seiner Leitfahigkeit x und 
seiner [relativen] Dielektrizitatskonstanten e. Die MAxweLlschen Gleichungen fiihren dann zu 
dem SchluB, da8 im Innern solcher Leiter kein elektrostatisches Feld existieren kann. Diese 
theoretische Voraussage widerspricht jedoch der Erfahrung: Neben dem vom Gitteraufbau 
fester Stoffe diktierten, rdumlich dreifach-periodischen Mikropotentialfeld treten im allgemeinen 
innerhalb ein und desselben Leiters merkliche Makropotentialdifferenzen auf, welche in den 
Kontaktspannungen zweier einander beriihrender Leiter manifest werden; die klassische Elektro- 
statik der leitfahigen Stoffe ist daher durch eine andere Auffassung zu ersetzen. 

Dem Aufbau einer diesen Anspriichen geniigenden Elektrostatik der Metalle wird die Fermi- 
SOMMERFELDsche Elektronentheorie zu Grunde gelegt: Das Gas der sogenannten freien Metall- 
elektronen gehorcht der spinmodifizierten Fermi-Statistik; es erweist sich im Grenzfalle ver- 
schwindender absoluter Temperatur als vollstandig entartet. Der dann bestehende Zusammen- 
hang zwischen der Konzentration jener Elektronen und dem elektrischen Makropotential wird 
kinetisch analysiert: Der statistische Gleichgewichtszustand des Elektronengases laBt sich als 
gegenseitige Kompensation eines vom Feld erzwungenen Konvektionsstromes und eines von Kon- 
zentrationsgradienten herriihrenden Diffusionsstromes genetisch verstehen; die je fiir diese 
Stréme maBgeblichen Werte der Beweglichkeit und des Diffusionskoeffizienten werden, nach 
Anpassung ihrer tiblichen Definitionen fiir merklich ideale Gase an den hier vorliegenden Fall des 
vollstandig entarteten Gases, explizit angegeben, und fiir ihren Quotienten wird eine einfache 
Beziehung ermittelt, welche als Ersatz der entsprechenden, auf den idealen Gaszustand be- 
schrankten EINSTEINschen Relation angesehen werden kann. 

Die Verbindung der FERMI-SOMMERFELDschen Elektronentheorie mit der PorissoNschen 
Raumladungsbilanz fiihrt auf eine nichtlineare, partielle Differentialgleichung zweiter Ordnung fiir 
das Makropotential, welche durch Einfiihrung , ,natiirlicher“‘ MaBeinheiten sowohl fiir das Poten- 
tial selbst wie auch fiir die Weglangen im Metall auf eine von den jeweiligen Daten des Metalles 
unabhangige Normalform gebracht wird; die genannten MaBeinheiten werden fiir das Beispiel 
des Kupfers zahlenmaBig berechnet. Die Differentialgleichung des Potentiales wird fiir parallel- 
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ebene Felder numerisch integriert, wobei drei Falle zu unterscheiden sind: Anreicherung des 
Elektrodenrandes mit Elektronen, Elektronenverarmung am Elektrodenrande, und vollige Ver- 
drangung der freien Elektronen aus einem Randgebiete, dessen Breite explizit angegeben wird; 
das Ergebnis wird durch Kurventafeln veranschaulicht. Im Lichte dieser Analyse erweist sich 
eine Revision der iitblichen Anschauungen iiber die Austrittsarbeit eines freien Elektrons aus 
seinem Muttermetall als unumganglich. Es wird vorgeschlagen, im Innern des Metalls begrifflich 
zwischen dem raumlichen Mittelwert des Potentiales y und dessen Erwartungswert <> zu unter- 
scheiden, welcher aus der Aufenthaltswahrscheinlichkeit der emissionsfahigen Elektronen an den 
verschiedenen Orten des Metalles resultiert ; die Differenz zwischen diesem Erwartungswert und 
dem radumlichen Mittelwert des Potentiales soll fiir die Austrittsarbeit verantwortlich sein. Dieser 
Vorschlag fiihrt auf eine neue Auffassung der sogenannten Feldemission: Mit zanehmender Feld- 
starke am Rande der Kathode steigert sich die ebendort herrschende Elektronenkonzentration, 
so daB das Einzelelektron an einem immer fiihlbarerem Platzmangel leidet; dieser geometrische 
Effekt zieht gemaB der Heisenbergschen Ungenauigkeitsrelationen die dynamische Zunahme 
seines Impulses nach sich, welche in einer entsprechenden Vergr6Berung der maximalen Null- 
punktsgeschwindigkeit jenes Elektrons kinematisch manifest wird. Daher existiert eine wohl- 
bestimmte, kritische Feldstarke, bei welcher die maximale Nullpunktsenergie der freien Elektro- 
nen am Kathodenrande deren Ubertritt in das angrenzende Vakuum gestattet. Diese ,,Durch- 
bruchsfeldstarke‘‘, bei welcher also die Feldemission eben einsetzt, stimmt in ihrer GréBenordnung 
mit jenen Feldstarken iiberein, bei welchen man bei der iiblichen Behandlung des gleichen Pro- 
blemes mittels des wellenmechanischen Tunneleffektes durch den Potentialberg der Grenze Me- 
tall-Vakuum auf merkliche Emissionsstromdichten gefiihrt wird. Umgekehrt 1aBt sich daher 
ein elektrostatisches Feld innerhalb eines planparallelen Plattensystemes nur dann verwirklichen, 
falls die Feldstaérke unterhalb der Durchbruchsfeldstaérke bleibt. Wird diese Bedingung ein- 
gehalten, so resultiert die Spannung eines parallelebenen Vakuum-Plattenkondensators aus der 
Vakuumspannung und den Eigenspannungen je der Kathode und der Anode; wahrend die 
Vakuumspannung stets proportional mit der Vakuumfeldstarke anwachst, besteht zwischen die- 
ser und jeder der beiden Elektrodenspannungen ein nichtlinearer Zusammenhang. DemgemaB 
bedarf der klassische Begriff der Kapazitat als einer spannungsunabhangigen KenngréBe des 
Kondensators einer Revision: Die Kapazitat erweist sich als Funktion der jeweils wirksamen 
Spannung. Ungeachtet der grundsatzlichen Bedeutung dieser Erkenntnis darf man sich praktisch 
meist mit dem Grenzwert jener Kapazitat bei sehr [infinitesimal] kleinen Absolutwerten der 
Spannung begniigen, dessen GréBe explizit angegeben wird. 


Franz Ollendorff, Haifa [Israel], Technion City El. Eng. Butlding. 
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Einleitung 


Ein fiir die Berechnung des Magnetisierungsbedarfes sehr wichtiges, bisher aber nur un- 
vollkommen geléstes Teilproblem des magnetischen Kreises ist die Beriicksichtigung des Feld- 
verlaufes im Luftspalt. Die herkGmmlichen Verfahren zur Bestimmung des Magnetisierungs- 
stromes basieren diesbeziiglich entweder auf ersten Naherungen und nicht immer ausreichend 
begriindeten Annahmen, oder es werden teilweise empirisch bestimmte Faktoren verwendet, 
deren Zuverlassigkeit in Grenzfallen nicht gesichert ist. Besonders unbefriedigend ist in diesem 
Zusammenhang auch die ungeniigende Beriicksichtigung des Jocheinflusses auf die Feldform. 
Offenbar entzieht sich die Behandlung dieses Problems besonders darum einer iibersichtlichen 
Lésung, weil der meist graphisch gegebene Verlauf der Magnetisierungskennlinie des Eisens 
zar Anwendung mathematisch wenig wirksamer Methoden verleitet. Da wenig Hoffnung 
besteht, ohne Anwendung des wirkungsvollen Hilfsmittels der Differentialgleichungen diese 
Aufgabe zu lésen, ist die erste Voraussetzung eine brauchbare analytische Nachbildung der 
Magnetisierungskennlinie des Eisens, wobei sicher ist, daB die Art der Nachbildung groBen 
EinfluB auf den rechentechnischen Aufwand und die mathematische Lésbarkeit iiberhaupt hat. 
Die vielfach geiibte Zuriickhaltung gegeniiber der Verwendung von Ersatzfunktionen fiir 
die Wechselstrom-Magnetisierungskurve wird meist durch Hinweis auf die Ungenauigkeit be- 
griindet. Zu beachten ist aber, da8 nur bei einer relativ kleinen Anzahl aller Vorausberechnun- 
gen die genauen magnetischen Eigenschaften des zu verwendenden Materials bekannt sind. 
Die innerhalb eines nicht kleinen Streubereiches méglichen Abweichungen von den mittleren 
Werten werden oft noch iibertroffen durch Ungenauigkeiten des sehr beengten Berechnungs- 
verfahrens, bedingt durch die graphisch gegebene sogenannte ,,genaue“‘ Magnetisierungskenn- 
linie. Ein besonderer Vorzug des analytischen Verfahrens ist neben der einfachen Anwendung 
seiner Ergebnisse der erheblich vergréBerte Informationsgehalt durch die allgemeinere und 
umfassendere Darstellungsméglichkeit der Ergebnisse. 


1. Ersatzfunktion fiir die Magnetisierungskurve 


Die Auswahl einer geeigneten Ersatzfunktion fiir die Wechselstrom-Magnetisierungskurve 
des Eisens erfolgt hier unter zwei Gesichtspunkten: Die Nachbildung soll eine ausreichende 
Genauigkeit vermitteln und im Hinblick auf die durchzufithrenden Rechnungen womdglich 
auf geschlossene Lésungen fihren. 

Bei empirisch gegebenen Magnetisierungskurven lassen sich die magnetischen Feldstarken 
durch mehrgliedrige Summen von verschiedenen Potenzen der magnetischen Induktion mit sehr 
hoher Genauigkeit darstellen [1]. 

H=khB+hk, BY +k, BY. (1) 


Diese Darstellung bietet im Vergleich zur inversen Formel B = B(H#) erhebliche rechentechni- 
sche Vorziige. In Ubereinstimmung mit dem iiblichen Vorgehen des Elektromaschinenbaues 
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soll unter B der Scheitelwert der Induktionsschwingungen und unter H der Effektivwert der _ 
magnetischen Feldstiarke verstanden werden. Eine befriedigende Nachbildung der Magneti- 
sierungskennlinie durch den Ansatz (1) erfordert mindestens ein Glied mit einer ziemlich hohen 
Potenz von B, wenn man wie in dem hier vorliegenden Fall auch den Sattigungsbereich der 
Kennlinie beschreiben will. Die obere Grenze fiir den Nachbildungsbereich, innerhalb dessen 
ausreichende Ubereinstimmung zwischen Magnetisierungskennlinie und Ersatzfunktion ver- 
langt wird, diirfte mit B = 18 bis 19 kG schon sehr vielen praktischen Belangen entsprechen. 
Bei Induktionen in der Nahe des Knicks der Magnetisierungskurve und darunter ergeben sich 
relativ kleine Feldstarken. Da die entsprechende magnetische Spannung im Vergleich zu der _ 
bei allen Asynchronmaschinen vorhandenen Luftspalterregung nur noch untergeordnete Be- 
deutung hat, brauchen die Genauigkeitsforderungen fiir den Bereich um 14 kG und darunter 
nicht zu hoch angesetzt werden. Fiir einen groBen Teil der Dynamoblechqualitaten scheint als 
einfachster und hinsichtlich seiner Genauigkeit annehmbarer Ansatz eine zweigliedrige Potenz- 
summe méglich, bestehend aus einem in B linearen Anteil und einer hohen Potenz in B (etwa 9 
bis 11). Die Gleichungen der folgenden Rechnung werden fiir einen solchen zweigliedrigen An- 
satz angeschrieben, obgleich sie auch fiir mehrgliedrige Potenzsummen (entsprechend Gl. (1)) 
erweitert werden kénnen. Der numerischen Lésung der Differentialgleichung wird fiir das 
zweite Glied als Beispiel die Potenz 9 zugrunde gelegt. Die Magnetisierungskonstanten sind fiir 
beide Glieder jedoch frei wahlbar. Sie gehen als Parameter in die Rechnung ein und erlauben 
damit von Fall zu Fall noch eine Anpassung an die gegebene Magnetisierungskennlinie. Sollten 
die Genauigkeitsforderungen trotz der frei wahlbaren Konstanten fiir auBergewéhnliche 
Magnetisierungskurven nicht befriedigen, so empfiehlt sich fiir die Potenzen ganze Zahlen zu 
verwenden und die Anpassung — wenn nétig — durch ein zusatzliches Glied zu verbessern. Die 
durchgefiihrte Rechnung basiert also auf folgendem Ansatz fiir die Magnetisierungskennlinie: 


He Boh Be (1a) 


2. Weitere Vorbereitungen und Voraussetzungen 


Der nichtlineare Zusammenhang zwischen magnetischer Induktion und Erregung fiihrt, 
wie schon angedeutet wurde, auf einen Feldverlauf im Luftspalt, der sich von dem des linearen 
magnetischen Kreises umsomehr unterscheidet, je weiter sich die Eisenmagnetisierung im 
Sattigungsgebiet befindet. Die Berechnung des Feldverlaufes laBt sich bei vorgegebenem 
Verlauf des Magnetisierungs-Strombelages mit guter Naherung durchfiihren, wenn man, der 
Vorschrift des elektromagnetischen Verkettungsgesetzes folgend, die fiir einen bestimmten 
Umlauf sich ergebenden magnetischen Spannungen und die elektrische Durchflutung ins Gleich- 
gewicht bringt. 


2.1 Magnetische Zahnspannung 


Bei der Berechnung der magnetischen Zahnspannung stellt sich zunachst die Frage, ob die 
Nachbildung (Gl. 1) der magnetischen Kennlinie auch zur Erfassung veranderlicher Querschnitte 
brauchbar ist. Bild 1 gibt eine bei elektrischen Maschinen gréBerer Leistung sehr haufig vor- 
kommende Zahnform wieder. Die Zahnbreite erweitert sich linear mit der normierten Koordi- 
nate ~ von der Breite % bis z,. Ist B,) die Zahninduktion an der Stelle w = 0, so ist mit der 
Zahnbreite z = Z (1 + ¢,+ u) an einer beliebigen Stelle w zwischen 0 und 1 die Zahninduktion 


Bai Bg 


°F an on u- (2) 

Oa 1 ist dabei die charakteristische Konstante der Zahnverjiingung. Die magnetische 
0 

Spannung, das Wegintegral tiber die magnetische Feldstarke, berechnet man fiir den Abschnitt a 

der Zahnlange /,, mit Gl. (1a) fiir die magnetische Erregung zu 


1 
du 
Via = bya | Bro hy [tae rok, | ool. (3) 
6 


Za 
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Die elementar ausfiihrbare Integration liefert fiir die beiden Integrale 


Ege s OE oy ieee aia 

Ca €q (m — 1) (1 + cq)" 
Die Gln. (3) und (3a) weisen auf die Méglichkeit hin, die geometrische Form des Zahnes durch 
modifizierte Magnetisierungskonstanten zu beriicksichtigen, etwa in der Art 


(3a) 


Cy,=h&2,, Ca =k, Lag» (3b) 
so daB man fiir die magnetische Zahnspannung des Abschnitts a schreiben kann: 
Vea = "24 (Bzo* Cig + B70" Cua) - (3c) 


Die beiden Faktoren Z, und Z, sind (Z,, fiir . = 9) im Bild 2 iiber der Verjiingungskennzahl c 
aufgetragen. 


jae ‘= 


we = 
20 0 0,2 0,4 0,6 0,8 4,0 1,2 44 1,5 
C— 
Bildi. Zur Berechnung der magn, Zahnspannung. Bild. Zur Berechnung der magn. Zahnspannung bei linear ver- 


Anderlichem Querschnitt, Z, linearer Ant., Z, nichtlinearer Ant. 


Besteht der Zahn aus zwei sich linear erweiternden Abschnitten a und 0, so kann man die 
magnetische Spannung berechnen aus 


Vogt Vay = Vz = lz (Bro: Gi + Bzo: C,) - (4) 


In dieser Gleichung wurde die gesamte Zahnlange 1, = 1,, + 1,, eingefithrt und die resultieren- 
den Magnetisierungskonstanten in folgender Form zusammengefaBt: 


lz . lzp lz Izy 
Cy — 7 Cia te 7 C1 , C,, om Coa es gay Cro 5 (4) 
Z Z Z Z 


C,, und C,,, sind analog zu C,, und C,,, zu bilden. Die gesamte Zahnspannung wird damit wieder 
auf die Form von Gl. (3c) gebracht. 

Es diirfte fiir die meisten Anwendungen ausreichend sein, die Zahnform durch aufeinander- 
folgende, linear veranderliche Querschnitte anzunahern, so daB man die Gl. (4) der magnetischen 
Zahnspannung fiir weitere Rechnungen zum Ausgangspunkt nehmen darf. Es ist dabei voraus- 
gesetzt worden, daB die Zahnentlastung durch parallelgeschaltete Nuten nicht beriicksichtigt 
zu werden braucht. Merkbare Fehler kénnen fiir Zahninduktionen iiber 18 kG entstehen, wenn 
das Verhaltnis Nut- zu Zahnbreite die Zahl 1 erreicht oder noch gréBer ist. Grundsatzlich ist 
die Beriicksichtigung der Zahnentlastung méglich; sie erfordert die Einfiihrung weiterer Para- 
meter. 


2.2 Magnetische Jochspannung 


Die Berechnung der Jochspannung erfordert die Erfassung eines zweidimensionalen magne- 
tischen Feldes. Solange Proportionalitat zwischen H und B besteht, erfolgt die Bestimmung 
des Jochfeldes durch Integration der Laplaceschen Differentialgleichung. In [2] wurde gezeigt, 
daB unter bestimmten Voraussetzungen fiir jeden Zusammenhang H = k 5” eine recht gute 
Naherung fiir das Jochfeld abgeleitet werden kann. Eine Voraussetzung bestand in der An- 


30 H. Wen: Analytische Behandlung des magnetischen Kreises von Asynchronmaschinen Bg tea 
Ean! 


nahme, daB die Permeabilitét des Joches durch die Umfangskomponente der Induktion be- 
stimmt sei. Wie sich zeigen laBt, ist in dem fiir die magnetische Jochspannung hauptsachlich 
maBgebenden Feldbereich diese Voraussetzung immer erfiillt. Eine weitere Voraussetzung ist 
die Kenntnis des értlichen Verlaufes der Induktion im Luftspalt. Damit kénnen die Induk- 
tionskomponenten des Jochfeldes als Ableitungen eines Vektor-Potentials berechnet werden, das 
sich in guter Naherung als Produkt zweier reiner Funktionen der Radial- bzw. der Umfangs- 
richtung ergibt. Die Verteilung der Umfangskomponente der Induktion iiber die Jochtiefe, 
abhangig von der Kriimmung des Joches, der Polzahl der Maschine, vom Verhaltnis Jochtiefe 
zu Polteilung und vom Kennlinienverlauf (Potenz m) kann berechnet werden. Bestimmt man 
die magnetische Jochspannung als Wegintegral der Feldstarke in Umfangsrichtung langs des 
Jochrandes am Nutengrund, so beschrankt sich das Feldproblem auf die Ermittlung der Um- 
fangskomponente der Jochinduktion an diesem Rand. Fiir sinusférmigen Induktionsverlauf in 
Umfangsrichtung berechnet man die auf den iiber die Jochtiefe gemittelten Wert der Induktion 
bezogene Randinduktion (y = y,) des AuBenjoches nach [2] zu: 


(ay 
=F lea es ol (5) 


Bir 


y, bedeutet den AuBenradius, 7, den Innenradius des jeweiligen Joches. AuBerdem ist 


—, i] in 2 
fae 4: / 1+ — 6) | (7) 
mit #, der Polpaarzahl. 


Vernachlassigt man die Kriimmung des Joches oder geht, was gleichbedeutend ist, zu un- 
endlich groBen Polzahlen iiber, so bestimmt man das entsprechende Induktionsverhaltnis fiir 


Y= 1 20 


Sa pa 
ha nas ctanh (2 ak (8) 


Fiir den Fall konstanter Permeabilitat (n = 1) wird aus Gl. (7) ¢4. = + p. Die Gln. (5) und (6) 
gehen damit in bekannte Beziehungen iiber. Auch Gl. (8) nimmt mit ” = 1 eine gelaufige 
Form an. 

Mit der Kenntnis der ungleichmaBigen Verteilung der Umfangskomponente der Jochinduk- 
tion laBt sich demnach die gesuchte Feldstarke am Nutengrund angeben zu 


Fiera ston (9) 


wo unter B; der tiber die Jochtiefe gemittelte Betrag der Umfangskomponente der Jochinduk- 
tion zu verstehen ist. In den Bildern 3 und 4 ist der Verlauf von s? fiir AuBen- und s% fiir 
Innenjoche jeweils fir = 1, 2,3 und s\ in Abhangigkeit von y/z; wiedergegeben. Bild 3 
zeigt die s-Werte fiir m = 1, Bild 4 fiir » = 9. Naherungsweise lassen sich die in den Gln. (5) 
und (6) angegebenen s-Werte auch zur Berechnung der Jochspannung fiir mehrgliedrige An- 
satze der Magnetisierungskurve nach Gl. (1) oder (1a) anwenden. Fiir den zweigliedrigen Aus- 
druck nach Gl. (1a) schreibt man die Umfangskomponente der Jochfeldstarke am Nutengrund 


. Be Bal acons } 
mit (10) 
IH, Jy = hy Sh - 


* Infolge eines Zahlenfehlers sind die im Bild 7 von [2] mit Gl. 8 fiir 7 = 11 und p > oo berechneten 
k,-Faktoren zu groB aufgetragen. 
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Die tiber die Jochtiefe veranderliche Verteilung der Feldkomponenten kann also naherungsweise 
in die Magnetisierungskonstanten der Ersatzfunktionen aufgenommen werden. Fiir die Berech- 
nung der magnetischen Jochspannung soll bei den folgenden Rechnungen von den fiir sinus- 


w= 
1,4 4 
2 Stander- 
joche 
3) 
co. 
| 
2 | Laufer- 
joche 
0,9 r 
4 
08 0,05 Oy 0/5 0,2 
2 BR 
Vite ot 
Bild 3. Zur Berechnung der magn. Jochspannung, linearer Anteil. 
0) 
p= 
4 
44 
| Stander- 
13 2 ¢joche 
3 
Ae 
| : 
Ur 44 
; 3 
40 = = = 2 \ Loufer- 
joche 
09 ; 
ost— ! i 
"0 0,05 0 0,5 0,2 
Ys IR 
she ae 


Bild 4. Zur Berechnung der magn. Jochspannung, nichtlinearer Anteil. 


formigen Induktionsverlauf im Joch bestimmten Faktoren s, und s, ausgegangen werden. Wie 
noch gezeigt wird, sind die Abweichungen des Induktionsverlaufes von der Sinusform im Joch 
wesentlich geringer als diejenigen der Luftspaltinduktion. Sie betragen fiir die hier zu Grunde 
gelegten Sattigungen im ungiinstigsten Fall etwa 5% bezogen auf die Grundwellenamplitude. 
Wie in [2] gezeigt wird, ist die Verteilung der Feldkomponenten tiber die Jochtiefe von der 
Feldverteilung in Umfangsrichtung abhangig. Der Einflu8 ist jedoch im Hinblick auf die er- 
wahnte geringe Felddeformation sehr klein und kann vernachlassigt werden. Hiernach berech- 
net man fiir geringe Abweichungen des Induktionsverlaufes von der Sinusform die Jochspan- 
nung nach der Gleichung 


t[2 


Tee ee a 
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wobei x die auf den Bohrungsdurchmesser bezogene Umfangskoordinate ist (Bild 5). Be- 
schrankt man B, ebenfalls auf Werte unter 18 bis 19 kG, so diirfen die durch parallel geschaltete 
Luftstrecken gebotenen Entlastungen unberiicksichtigt bleiben. Bei an das Lauferjoch an- 
grenzender Welle soll die Entlastung in erster Naherung durch die Wahl einer entsprechend 
vergroBerten Lauferjochdicke erfaBt werden. 

Nachdem das zweidimensionale Jochproblem niherungsweise auf eine eindimensionale 
Aufgabe, namlich die Integration der Gl. (11) zuriickgefiihrt ist, steht der Ermittlung 
des Feldverlaufs in Umfangsrichtung nichts mehr im Wege. 


3. Feldverlauf und gesamter Magnetisierungsbedarf 
3.1 Die Differentialgleichung des Feldverlaufes 


Bild 5 gibt zur Orientierung das Schema eines magnetischen Kreises wieder. Auf der linken 
Bildhafte ist eine mégliche Form von Zahnen und Nuten der Wirklichkeit entsprechend an- 
gedeutet. Mit Hilfe der im Abschnitt 2.1 angegebenen Gleichungen wird es méglich, hierfir 
Ersatzzahne konstanter Breite (rechte Bildhalfte) anzunehmen, wenn die Magnetisierungs- 
konstanten in der Nachbildung der magnetischen Kennlinie entsprechend geandert werden. 
Ausgehend von Gl. (4) schreibt man die 
magnetischen Spannungen fiir Stander- 
(Index S) und Lauferzahne (Index R) 
entsprechend den Gln. (12a), wobei die 
magnetische Spannung und die Zahn- 
induktion als stetige Funktionen von % 
zu deuten sind. Mit den im Abschnitt 2.2 
angegebenen Abkiirzungen, mit Gl. (11) 
und unter den genannten Voraussetzun- 
gen berechnet man die magnetische 
Spannung, innerhalb des Polbogens o 
bis x fiir Stander und Lauferjoch nach 
den Gln. (12b). 


Vis =lzs(Cis + Bzs + Cas: BZs), 
Ver=lzr(Cir: Ber t+ Curr: Ber) 


| (12a) 


| Uiorclae 28 fee is Palas Bis] dx, 


te) 


| lp SS i [Use Bie fae Bj] dx, 


Bild 5. Magnetischer Kreis, Koordinatensystem und Bezeichnungen. 


ee 


(12b) 


t;s bzw. tr sind die Polteilungen langs des Integrationsweges fiir die Jochspannung (siehe 
Bild 5), wahrend 7 die Polteilung langs der Bohrung ist. Das negative Vorzeichen folgt aus der 
Wahl des Nullpunktes fiir x im Maximum des Luftspaltfeldes und der Forderung nach einem 
positiven Wert fiir V;; und V;z in den Gln. (12b). AuBerdem ben6tigt man fiir die durchzu- 
fiihrenden Reenanseen die Beziconune zwischen Zahninduktion und Luftspaltinduktion, sowie 
den Zusammenhang zwischen Jochinduktion und Luftspaltinduktion. Diese vier Gleichungen 
sind gegeben in folgender Form: 


(13a) 


Te | By: dx, Bir = Bys. =. (13 b) 
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In den Gln. (13a) bedeutet ¢ die Nutteilung und / einen Faktor, der die Eisenfiillung be- 
riicksichtigt. Zur Erfassung von Streufliissen im belasteten Zustand sind die z- und y-Werte 
entsprechend zu reduzieren. Denkt man sich in der ersten Gl. (13b) die Luftspaltinduktion 
B, zerlegt in Grund- und Oberwellen, so wird verstandlich, daB der Verlauf der Jochinduktion 
einen erheblich geringeren Oberwellengehalt aufweist, als der Verlauf der Luftspaltinduktion. 
Einer Oberwelle der »-ten Ordnung im Luftspaltfeld mit dem Betrag 1 entspricht in der 
Umfangskomponente der Jochinduktion eine Oberwelle gleicher Ordnungszahl vom Betrag 1:». 
SchlieBlich verbleibt noch die Berechnung der magnetischen Spannung des Luftspaltfeldes. 
Mit 6, dem die Nutung beriicksichtigenden aquivalenten Luftspalt konstanter Breite, be- 
rechnet z. B. durch Anwendung der Carterschen Luftspaltfaktoren, wird 
V,=<8,, (14) 
Mo 


wobei fy die Permeabilitat des Vakuums ist. 
Gibt man wie tiblich den Magnetisierungsstrombelag sinusférmig vor, setzt also 


a(x) = — A,- sin (= *) ; (15) 


so folgt aus dem elektromagnetischen Verkettungsgesetz, das die Gleichheit der magnetischen 
Spannungen und der elektrischen Durchflutung fordert 

eee te Ve ie ee Ay-sin(2 x) ax (16) 

T 

Eine einmalige Differentiation nach x, sowie die Verwendung der Beziehungen (13) und (14) 
zwischen den magnetischen Spannungen und den Induktionen fiithrt nach einigen Umformungen 
auf eine gewohnliche nichtlineare Differentialgleichung zweiter Ordnung fiir die Feldverteilung. 
Es empfiehlt sich als Veranderliche die Jochinduktion B;, einzufiihren, so daB folgt: 


” (7 rr ” % A ae 
Vof- Bjs + aN (ys ies (Bjs) : Bs as Bjs is SS ap Sin (= x). (17) 


Die gestrichenen Gré8en bedeuten in Gl. (17) Ableitungen nach der Umfangskoordinate »%. 
AuBerdem wurden folgende Abkiirzungen eingefiihrt: 


Br eee oe en | 
vine Ta tes Cns(7) + ler Cur (32 


) M,= 2 e Jes "i (22 y: Fu (17a) 


, 


as US ES 
M, =F [Elst 2 he 
R 


D ce Tt D T YR 
ey fa Ga os Pe Oe 
= — + bys Cis —— + ereir=-:; 
Ho I 25 S2R 


3.2 Loésung der Differentialgleichung (17) 


3.2.1 Vernachlassigung der Jochspannungen 


Fiir den Sonderfall vernachlassigbar kleiner magnetischer Jochspannungen, also fiir y; 00 
und y, 00, wird aus Gl. (17) und der Substitution. 


is fY¥s E (18) 
, eae van ig Ihe 
Bi(1 + aN - By) = —Asin (= x). (19) 
Die Dgl. (19) laBt sich integrieren und man erhalt fiir die Luftspaltinduktion die Losung 
B, +N Br = Bycosa, (20) 
res 
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wobei zur Abkiirzung By = oe und « = ~ xeingefiihrt wurden. Die Integrationskonstante 
It tT 
ergibt sich zu 0; sie wird aus der Bedingung B, = 0 fiir ¥ = = bestimmt. 


Die Gl. (20) beschreibt also den Feldverlauf im Luftspalt bei vernachlassigbar kleinem Joch- 
einflu8. Man kann in dieser Gleichung Bo, das die Einheit einer magnetischen Induktion hat, 
als Ma fiir die GréBe des Magnetisierungsbedarfes ansehen. 


OA 0,2 04 06 40 48 S,=20 


4,0 
0,9 = = 
08 ss See as 
SS 
0,7 ~ < — 
06 
ae Se a eee 
00 40 70 700 
<efeannn 


Bild 6. Graphische Ermittlung von #, als Funktion von S,. 


AG oe iz == 2 
; 4 yo 
2.0 i 1 
S 4 
48 at i 74-005 
SS 1 
1,6 > | 
RS 
=0 
44 4 : 4 
42° 
41,0 


Os ek? Oh Ol — 0s 24  eek  atee nn eee eon 


4 


Bild 7. Zur Berechnung des Magnetisierungsbedarfes in Abhangigkeit von Joch- und Zahnbelastung, 
verschiedene Naherungen mit und ohne Beriicksichtigung der 3. Oberwelle der Jochinduktion. 


Sieht man als Ziel der Aufgabe nicht nur die Bestimmung des Feldverlaufes im Luftspalt, 
sondern auch die Ermittlung des Magnetisierungsstrombelages, so ist By als unbekannte GréBe 
anzusehen. Als bekannt darf man zu ihrer Ermittlung die Amplitude der Grundwelle der 
Luftspaltinduktion 6,, vorlegen. 

Zur Bestimmung von A, kann folgender Weg eingeschlagen werden: Dividiert man Gl. (20) 
durch By und bezeichnet B = B,/Bo, fiihrt auBerdem den Zahnparameter S = Bj—*- N ein, 
wodurch sie iibergeht in 


pas Bo = cose, (21) 


so 1aBt sich der Kurvenverlauf B(x) fiir S = const und damit auch nach Ausfiihrung einer 
FourtgRanalyse fp, = B,,/By in Abhangigkeit von S angeben. f, = B,(S) ist im Bild 6 ein- 
gezeichnet. Schreibt man den Zahnparameter S mit Hilfe von B, um in 


S =p) = yy Une = See ee (22) 


tragt auBerdem S fiir verschiedene Parameter S, im gleichen Bild 6 auf, so sind die Schnitt- 
punkte der Kurven Wertepaare des gesuchten Zusammenhanges f, = f,(S,). Anstelle von p, 
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0 


13% 1 : i : 5 
2 ce Geaptte Abhangigkeit von S, aufgetragen. Der Magneti- 
1 L1 


sierungsstrombelag berechnet sich dann nach Gl. (23) zu 


ist im Bild 7 (Parameter y, = 0) 


A B 
~=—* (1+4,;), 
mit Ay, ae und —— Izs+Cys ath + lz p Cee ; 
8 Mo 0) fz St ZR 


A, ,. bezeichnet jenen Strombelag, der bei unendlich groBer Leitfahigkeit des Eisens aufgebracht 
werden mu8. d; erfaBt das Verhaltnis des linearen Anteils der Zahnspannungen zum Luft- 
spaltanteil der magnetischen Spannung und ist klein gegen 1. 


3.2.2 Sinusformiger Verlauf der Jochinduktion 


Die Integration der Dgl. (17) ist nicht mehr in geschlossener Form méglich, wenn die magne- 
tische Jochspannung beriicksichtigt werden soll; man ist also zu ihrer Lésung auf die Zuhilfe- 
nahme numerischer Methoden angewiesen. Um einen Uberblick iiber einen zweckmaBig er- 
scheinenden Gang der Lésung zu geben, werden als erste Naherung Feldverlauf und Magnetisie- 
rungsbedarf berechnet, wie sie sich bei Beriicksichtigung der Grundwelle der Umfangskompo- 
nente der Jochinduktion ergeben. In Gl. (17) setzt man also in den mit M multiplizierten Glie- 
dern den angenaherten Verlauf der Jochinduktion 


B;s = Bys,-sine (24) 
als bekannt voraus. So wird dann 


Bis = Bes,- = [126 sin « — 84 sin 3a + 36sin5a—gsin 7a + singa]. (24a) 


Schreibt man die beiden Jochglieder auf die rechte Seite der Dgl. (17), fiihrt wieder die Luft- 
spaltinduktion B, ein und integriert beide Seiten, so entsteht 


BON BY = (ies M, Bjs, — M,, B}s,- 0,492] cosa + Moe B?s,+0,1094+ B(a), 
ise . : i (25) 
mit W(x) = cos 3 a — 0,257 cos 5 « + 0,0459 Cos 7 4 — 0,00397 Cos ga. 
Kiirzt man weiter ab 
; A 
By = — eee! Bis, — M, Bi s,-04492, | 
SN Se an | ae a Se a (25) 
rae B M,,— Bis, 0,1094 = 1 aie vi° Bs | 
Set ogee tea Ge 
so kann man Gl. (24) auf die fiir numerische Rechnungen handlichere Form bringen: 
B+ S’- B® =cosa+y- O(a). (25b) 


Wahlt man zunichst y als Jochparameter (y = o entspricht der bereits bekannten Lésung), so 
kann nach dem oben schon beschriebenen Lésungsverfahren vorgegangen werden. Das Er- 
Bri 
Bo 
abmessungen und der Grundwelleninduktion B;, abhangt, wahrend y noch die zu suchende 
GroBe B, enthalt. In der zweiten Gl. (25a) ist bereits angedeutet, wie der Zusammenhang zu 
dem von der Lésung unabhangigen Parameter y, ist. Man kann danach, wenn f, bekannt ist, 
jedem y ein y,, das kennzeichnend fiir den JocheinfluB ist, zuordnen. 

Im Bild 7 sind einige das Ergebnis der Rechnung zusammenfassende Kurven fiir 1 = 9 
und einen gewissen Bereich der Hilfsveranderlichen S, und y, angegeben. Die aufgetragenen 


Bie 


= B,(S;, y), worin S, wie oben nur noch von den Maschinen- 


gebnis ist eine Funktion pf; = 
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Zahlenwerte B,/B,, sind durch die erste Gl. (25a) definiert. Die Doppelgl. (26) gibt den Zu- 
sammenhang mit dem Strombelag an (linke Seite) und weist, durch Einfiihrung von 7, auf 
den Berechnungsgang hin (rechte Seite). 
B ath 72h B 
eee & M, Biss) = 52 +45%:- (26) 


Bry, Bin 


B,/Br, kann also berechnet werden, sobald fp, = B, ,/By. graphisch in Abhangigkeit von S, 
und y, gegeben ist. Liegt B,/B,, vor, so kann aus Gl. (26) der resultierende Strombelag, be- 
zogen auf A,., entsprechend Gl. (23) ermittelt werden zu 


lo 

Tt Bis 

= LES 2 
Aico Bry 14 Bry ( 7) 

Der bezogene Magnetisierungsstrombelag dieser Gleichung wird maBgeblich bestimmt durch 

die beiden Hilfsveranderlichen S, und y, die den nichtlinearen Anteil der magnetischen Zahn- 

und Jochspannung erfassen. Mit Hilfe der Gln. (22) und (17a) kann man schreiben: 


Cree Cas: Basilzs+ Car Bzrilzr (28) 
0: Bry 
Mo 

Unter Bz;, und Bzp, sind die der Amplitude der Grundwelle der Luftspaltinduktion B,, 

entsprechenden Induktionen fiir Stander- und Lauferzahne (an deren engster Stelle) zu ver- 

stehen. Mit der erwaihnten Bedeutung der Magnetisierungskonstanten C ist S, durch Gl. (28) 

als Verhaltnis der nichtlinearen Anteile der magnetischen Zahnspannungen zum linearen 

Anteil von Luftspalt und Zahnspannung, zu berechnen fiir die Grundwellenamplitude der 

Induktion, kenntlich gemacht. Dem Parameter y, kommt nach den Gln. (25a) und (17a) eine 
ahnliche Bedeutung zur Erfassung des Jocheinflusses zu. Es ist 


+ Cys Bzsilzs + Cir Bzri!zr 


9 9 
Jns*Bysitys t+ Jnr Bpritir 


71 = 0,0348 ~ (29) 


A + (4s Bzsilzs + Cir Bzri!zr 

0 

y, setzt einen Teil der nichtlinearen magnetischen Jochspannungen ins Verhaltnis zu den line- 
aren Anteilen von Luftspalt und Zahnspannungen. Alle Anteile sind wieder fiir die Grund- 
wellenamplituden der Luftspaltinduktion zu berechnen. Es ist also z. B. 


Bis, = Br. (29a) 


Wiirde man, wie meist iiblich, mit Hilfe von k,-Faktoren rechnen, die definiert sind nach der 
Gleichung 


qt; 


so lieBe sich auch fiir den nichtlinearen Anteil der magnetischen Jochspannung ein entsprechen- 
der k,-Wert ermitteln zu 


n/2 


So: i sin? & da = 0;2595). (31) 


19) 


Im Faktor s} wiirde entsprechend Bild 4 der Einflu8 von Kriimmung, Jochtiefe und Polzahl 
zum Ausdruck kommen. Fiir y, k6nnte man dann schreiben: 


hykigs a Bossy +- Ry hy9R = Biri 
= 0,269 ——— = = 
V1 , ) 0 Bri aN 7 (32) 
aan ce Cs Bzslzgs+QyrBerizr 
0 
Der vor dem Bruchstrich stehende Zahlenwert rechtfertigt die Bemerkung, daB der Zahler 
von y nur einen Teil der nichtlinearen Spannungsanteile darstellt. Eine auf volle Erfassung 
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der magnetischen Spannungsanteile abhebende (MaBstabs-)Anderung fiir y, bringt jedoch keine 
weiteren Vorteile, da es auch dann nur eine fiktive GréBe ist. Joch- und Zahnspannungen k6n- 
nen nicht durch y, oder S, allein erfaBt, sondern sie miissen jeweils durch S, und y, gemein- 
sam beschrieben werden. Es wird deshalb weiterhin y, nach Gl. (29) bestimmt. 

Durch die Vernachlassigung vor allem der 3. Oberwelle in der Umfangskomponente der 
Jochinduktion entsteht ein gewisser Fehler. Er ist bei kleinen S,-Werten und gréBeren Joch- 
beanspruchungen so gerichtet, daB das vorher beschriebene Verfahren etwas zu groBen Magne- 
tisierungsbedarf angibt. Bei groBen S,-Werten berechnet man etwas zu kleine magnetische 
Spannungen. 


3.2.3 Feldverlauf mit Beriicksichtigung der dritten Oberwelle in der Jochinduktion 


Die folgende (zweite) Naherung fiir die Verteilung des Luftspaltfeldes beriicksichtigt auch 
die dritte Oberwelle der Umfangskomponente der Jochinduktion. Eine Heranziehung weiterer 
Oberwellen ist einerseits durch den rasch zunehmenden rechnerischen Aufwand begrenzt, 
andererseits eriibrigt sich die Erfassung von Oberwellen héherer Ordnungszahlen aus den schon 
angedeuteten Griinden. Wie die Analyse der Feldkurve bei Beriicksichtigung der Grundwelle 
der Jochinduktion zeigt, erreicht das Verhaltnis der Luftspaltinduktionen B, ,/B,, innerhalb 
des untersuchten Bereiches von S, und y,; nur an einer Stelle den Betrag von 5°/). Die 
fiinfte Oberwelle tritt also in der Umfangskomponente der Jochinduktion im ungiinstigsten 
Fall mit 1% der Grundwellenamplitude in Erscheinung, so daB eine Vernachlassigung bereits 
als tragbar angesehen werden kann. Dies gilt in noch starkerem MaBe fiir Oberwellen héherer 
Ordnungszahlen. 

Fir die anschlieBenden Rechnungen sei das Verhaltnis der Amplitude der Jochinduktion 
mit 

Bis3 = 1 Brg =} 


Bjs1 3) Bina (33) 


bezeichnet. ¢ wurde durch FourIERanalyse aus dem Feldverlauf bei ausschlieBlicher Beriick- _ 
sichtigung der Grundwelle der Jochinduktion gewonnen und bildet den Ausgangswert fiir 
weitere Naherungen. In die Dgl. (17) setzt man nun fiir die Anteile der Jochinduktion 

B;s = B;s, (sina —#sin 34). (34) 


Fiir B?, erhalt man einen formelmaBig etwas umstandlichen Ausdruck, der hier nur abgekiirzt 
wiedergegeben wird, in der Form 


(fs) T, sino — T,sin30+T7,;sin5a—T,sin7«4 + Tysinga. (35) 
jS1 


-01 008 006 -0,04 0,02 0 66002 ©6004 006 008 of 


Bild 8. Zur Berechnung von BY mit Beriicksichtigung der dritten Oberwelle. 


Die Koeffizienten 7, bis T, berechnet man aus einer Potenzreihe, die Glieder bis # enthalt. 
Fiir den Bereich —0,1 <¢ < 0,1 sind die Koeffizienten im Bild 8 wiedergegeben. Setzt man 
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die Ausdriicke (34) und (35) in die Dgl. (17) ein, so 1aBt sich diese wieder integrieren und man 
erhalt folgendes Ergebnis: 

dy, 
3 


i 


Bo By ee P, (a, t) (36) 


= My Bg, 2M Bie, T,] 005 a + M,— Bis, 


Hierbei wurde das Glied M,¢ B,s, als sehr klein gegeniiber den anderen Anteilen vernach- 
lassigt. Kiirzt man analog der Gl. (25a) wieder ab, 


, tee f 
Boe |S ae ea eee 1), 
az} D 
ee a (36a) 
Mie Swag ‘ip 
lie ys 3 Bo: ela Bi 00,3282 ” 


so kann Gl. (36) auf die der Gl. (25) entsprechendeForm gebracht werden: 

Bp + S’ p? = cosa + y,- O(a, 7) . (37) 
Ihre Auswertung kann ahnlich wie bereits beschrieben erfolgen. Im Bild 7 ist das Verhaltnis 
B,/Br, wie fiir Gl. (26) in Abhangigkeit von S, und , die identisch mit den dort angegebenen 
Gr6éBen sind, aufgetragen (2). Es wurde berechnet iiber die Beziehung 


af, 
BP; We es (38) 
Bry Bry 0,1094 


Die Berechnung des Strombelages erfolgt wieder nach Gl. (27). 


3.2.4 Zur Konvergenz des Verfahrens der schrittweisen Naherung 


Die im Bild 7 sichtbaren Unterschiede im Magnetisierungsbedarf, verursacht durch die in 
erster Naherung beriicksichtigte dritte Oberwelle der Jochinduktion machen sich besonders 
bei hohen Jochbeanspruchungen bemerkbar. Man berechnet bei ausschlieBlicher Beriicksich- 
tigung der Grundwelle der Jochinduktion und hohen Zahnbeanspruchungen zu groBe, bei klei- 
nen Zahnbeanspruchungen zu kleine magnetische Erregungen. Die schrittweise Naherung zur 
Bestimmung der dritten Oberwelle ist im Bild 9 verdeutlicht. Mit 1 sind jene ¢-Kurven ge- 
kennzeichnet, die sich aus der ersten Naherung bei Vernachlassigung der Oberwellen der Joch- 
induktion (fiir ¢= 0) ermitteln lassen. Sie sind aufgetragen iiber y, mit dem Parameter der 
Zahnbeanspruchung S = const. Die aus der ersten Naherung berechneten ¢-Werte dienen 
zur Bestimmung der zweiten Naherung fiir Feldverlauf und Magnetisierungsbedarf. Daraus 
resultieren die bezogenen Amplituden der dritten Oberwelle, die mit 2 gekennzeichnet sind 
(2. Naherung). Sie hegen zwischen dem Wert ¢ = 0, also der nullten Naherung und den 
t-Kurven der ersten Naherung. 

Zieht man das Ergebnis der zweiten Naherung fiir ¢ zur Berechnung der dritten heran, so 
findet man die entsprechenden (mit 3 gekennzeichneten) ¢-Werte wieder zwischen den Kurven 
der beiden vorangegangenen Naherungen. Die Konvergenz der schrittweisen Naherung wird 
hierdurch deutlich. 

Ahnlich wie im Bild 9 die beiden vorhergehenden Lésungen die folgende eingrenzen, ist 
dies auch bei der Darstellung des Magnetisierungsbedarfes B)/B,;, im Bild 7; erste und zweite 
Naherung grenzen die dritte und alle folgenden ein. Man beschleunigt das Konvergieren des 
Verfahrens der schrittweisen Naherung in dem man z. B. zur Berechnung der dritten Naherung 
t-Werte vorlegt, die (zwischen den Werten der ersten und zweiten Naherung liegen und) den 
zu erwartenden endgiiltigen #-Werten méglichst gut entsprechen. Ob nach einem weiteren 
Schritt nochmals Anderungen des Magnetisierungsbedarfes und des Feldverlaufes zu erwarten 
sind, zeigt ein Vergleich der Ergebnisse zwischen n-ter und (m + 1)-ter Lésung. Die so als 
dritte Naherung berechnete Loésung fiir den Magnetisierungsbedarf ist fiir eine etwas feinere 
Unterteilung des Par ameters y,1m Bild 10 aufgetragen, ihre Lage im Vergleich zur ersten und 
zweiten Naherung ist im Bild 7 (gekennzeichnet durch 3) angedeutet. Bild 9 zeigt ebenfalls die 
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zu dieser Losung gehorigen ¢-Werte (3). Die bei der Fortfithrung des Verfahrens zu erwarten- 
den endgiiltigen Lésungen unterscheiden sich von den Werten der angegebenen dritten 
Naherung an der ungiinstigsten Stelle (Sy = 2,0, y, = 0,15) mit Sicherheit um weniger als 
2% und an allen anderen Stellen um noch kleinere Betrage, so daB die aufgetragenen GroBen 
als gute Naherung fiir die endgiiltige Losung angesehen werden kénnen. 


—>t 


~0,025 


~0,05 4 


> 
0,05 0" 015 
etc 


Bild 9. Anteil der dritten Oberwelle in der Jochinduktion, verschiedene Naherungen. 


(x) 
a4 


22 


iar te 06 = Suet 42. F146 +18 ~~ ~2,0 


yi 


Bild 10. Zur Berechnung des Magnetisierungsbedarfes in Abhangigkzit von Joch- und Zahnbelastung, 
S, und y, werden mit Grundwellen-Induktionen bestimmt. 


Da mit den beiden Hilfsveranderlichen S$, und y, das nichtlineare Verhalten des magneti- 
schen Kreises vollstandig beschrieben wird, kann man auch die sogenannte ,,Abflachung“ oder 
allgemeiner, die Deformation des Luftspaltfeldes (an der Stelle a = 0) in Abhangigkeit von 
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diesen beiden Werten beschreiben. Im Bild 11 ist das Verhaltnis der maximalen Luftspalt- 
induktion B, max bezogen auf den Scheitelwert der Grundwelle B,, iiber der Zahnsattigungs- 
kennzahl S, mit dem Parameter y, = const fiir die Jochsattigung aufgetragen. Fir y, > 0 
ergeben sich im Bereich kleiner $,-Werte Zuspitzungen der Feldkurve gegeniiber dem sinus- 
formigen Verlauf, also das Gegenteil der Abflachung. Die Kurven bestatigen, daB sich Zahn- 
und Jocheinflu8 im Hinblick auf die Felddeformation in gegenlaufigem Sinne auswirken. 


4 naan Gee 
ee 
4,0 
Ss ese Y= 
0/15 
0,9 0; 
0,05 
0 
0,8 eae 
0 0,2 0/4 0,6 08 4,0 12 44 1,6 1,8 2,0 


4 


Bild 11. Maximalwert zu Scheitelwert der Grundwelle der Luftspaltinduktion in Abhangigkeit von Joch- und Zahnbelastung, 


Zusammenfassung 


Bildet man die empirisch gegebene Magnetisierungskurve H = H(B) als Summe mehrerer 
verschiedener Potenzen von B nach, so 1aBt sich ihr charakteristischer Verlauf bereits durch 
wenige Glieder wiedergeben. Die Berechnung von Feldverlauf und Magnetisierungsbedarf 
wird hier auf der Grundlage einer durch eine Potenzsumme dargestellten Magnetisierungs- 
kennlinie entwickelt. Es zeigt sich, daB das nichtlineare Verhalten des magnetischen Kreises 
durch zwei Parameter beschrieben werden kann, die aus den nichtlinearen Anteilen von Zahn- 
und Jochspannung, bestimmt fiir die der Grundwellenamplitude entsprechenden Induk- 
tionen, berechnet werden. 

Fiir die durchgefiihrten numerischen Rechnungen zur Ermittlung von Feldverlauf und 
Magnetisierungsbedarf wurde die Magnetisierungskurve durch ein lineares Glied und ein 
Glied mit einer neunten Potenz in B nachgebildet. Die Wahl der Konstanten dieser Glieder 
ist beliebig und kann dem jeweils vorliegenden Fall angepaBt werden. 
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Die Steuerkennlinien eines Stromrichters 
mit antiparallel geschalteten Ventilen bei Wechselstrombetrieb 


Von 
Kurt MULLER-LUBECK, Mannheim 


Mit 9 Textabbildungen 


(Eingegangen am 3. Juni 1960) 


Ubersicht: Es werden die Steuerkennlinien eines Stromrichters mit antiparallel geschalteten Gasent- 
ladungsgefaBen bei Wechselstrombetrieb mit einer aus einem Ohmschen Widerstand und einer Induktivitat 
bestehenden Belastung beschrieben. 


1. Aufgabestellung 


Sehr bald nach der vor mehr als 30 Jahren erfolgten Entwicklung gittergesteuerter Gas- 
entladungsventile, also der Einrichtung von Mehranoden-GleichrichtergefaBen oder Einanoden- 
Quecksilbersampf-Gliihkathodenréhren mit Steuergittern, wurde die Anwendung derartiger 
Stromrichtergefa8e nicht nur fiir den Gleichstrom-, sondern auch fiir den Wechselstrombetrieb 
in Betracht gezogen. Wahrend der Gleichstrombetrieb bekanntlich dem Gleichrichter, dem 
Wechselrichter und den damit gebildeten Umkehrschaltungen vorbehalten ist, liegt ein 
Wechselstrombetrieb bei den Umrichtern und bei den Schaltungen fiir eine einfache Wechsel- 
stromtibertragung vor. 


Der letztgenannte Wechselstrombetrieb von Stromrichtern in verschiedenen Schaltungen 
ist schon frithzeitig fiir die Steuerung von Punkt- und NahtschweiBmaschinen [1], [2], [3] und 
zur Helligkeitssteuerung von Beleuchtungsanlagen angewendet worden. In neuerer Zeit 
sind derartige Stromrichter fiir Wechselstrombetrieb mit einanodigen GefaéBen fiir hohe 
Stromstarke und hohes Sperrvermégen als Wechselstromschalter und auch zur Lésung von 
Erwarmungsproblemen in der Hiittenindustrie, beispielsweise fiir die sog. unmittelbare Wider- 
standserwarmung von Kniippeln bei direktem Stromdurchgang [4], interessant geworden. 


Die Steuerung der Stromrichter bei Wechselstrombetrieb ist, wie sich zeigt, an den Liick- 
betrieb der Ventile gebunden. Uber diesen Liickbetrieb liegen rechnerische Unterlagen nur 
bei Gleichrichteranwendungen vor [5]. Deshalb soll nachfolgend die Berechnung der Steuer- 
kennlinien dieser Stromrichter im Wechselstrombetrieb mit antiparallel geschalteten Gas- 
entladungsréhren nachgeholt werden. Die Untersuchung geht von einer einfachen Grund- 
schaltung eines Wechselstromkreises mit antiparallel geschalteten gesteuerten Ventilen und 
der Annahme einer aus einem Ohmschen Widerstand und einer Induktivitat bestehenden 
Belastung aus. 


2. Betriebsweise, Steuerbereich und Stromflu@dauer 
Die Grundschaltung zeigt Bild 1. Die Wechselspannung u,, habe den zeitlichen Verlauf 
u, =V2U,sinx, «=o, (1) 


die Belastung bestehe aus dem Widerstand R und der Reaktanz X = w L. Die Ventile sind 
mit V, und V, bezeichnet. Sind beide Ventile standig stromdurchlassig, so besteht der normale 
Wechselstrombetrieb, fiir den der Momentanwert des stationaren Wechselstromes 


D 


w RR? i x2 Z 


a et 
__}? Uae ie (x —o) = Ee e050 sin (x —o), | 


mit (2) 
tt | 
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betragt. Der Wechselstrom ist wie die Spannung sinusférmig, hat den Effektivwert 


a4 Uy _ Uy 
dae ees F608 @ (3) 


und die Phasenverschiebung @ gegen die Spannung. 


Es soll untersucht werden, welche Verhiltnisse bei Gittersteuerung der Ventile eintreten’ 
Sie erfolge durch Gitterspannungen, die die Ventile im Ziindaugenblick x = a bzw.% =a+a2 
stromdurchlassig machen. Die dazu erforderlichen Gitterspannungsimpulse werden aus dem 
Gittersteuergerat St geliefert. Da der natiirliche Ziindeinsatz ohne Gittersteuerung bei standig 
stromdurchlassigen Ventilen bei x = oe 
bzw. *=o+ a liegt, treten offenbar 
verschiedene Strombedingungen ein, ie 
nachdem « < @ oder a > @ ist. 


ly 
<+—_ > 


ican 


R X-ol 


Bild 1. Wechselstromkreis mit antiparallel geschalteten Strom- 
richtergefaBen mit einer aus einem Ohmschen Widerstand 
und einer Induktivitat bestehenden Belastung. Bild 3. Strom- und Spannungsverlauf bei « > @. 


a) Zindwinkel « <oe 


Dieser Fall ist in Bild 2 veranschaulicht. Die Freigabe des Ventiles V, erfolgt, waihrend 
das Ventil V, noch stromfiihrend ist. Die Ventilablésung erfolgt ohne besonderen Kom- 
mutierungsvorgang im Nulldurchgang des Stromes. Die StromfluBdauer 4 jedes Ventiles 
betragt A = a. Der Effektivwert des Stromes hat den in Gl. (3) angegebenen Wert unabhangig 
von «. Die Anodensperrbeanspruchung der Ventile ist Null. 


b) Zindwinkela>o 


Bei einem Ziindwinkel « > @ tritt eine Verspadtung des Ziindeinsatzes gegeniiber x = 0 
und infolgedessen eine verkiirzte StromfluBdauer A und somit liickender Strom ein. Es findet 
damit eine Herabsteuerung des effektiven Wechselstromes statt. 


Zur Berechnung des Wechselstromes hat man von der Spannungsgleichung 
nee 


Uy =X a ae Aaa 


auszugehen, zu der die Anfangsbedingung 7, =o fiir x =a hinzutritt. Dies fiihrt nach 
_ Integration der Differentialgleichung in Verbindung mit der Anfangsbedingung zu der Strom- 
gleichung 


# cos @ (sin (x — e) — sin (« — @) e—7 2) , (4) 


die fiir x = « bisx =a +2 gilt. 
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Bei x = « +A erfolgt der Nulldurchgang des Stromes und damit das Verléschen des zu- 
gehérigen Ventiles. Der Wert 4 der StromfluBdauer ergibt sich somit aus der Bedingung 


sin (a —o +A) — sin (a — 9) ¢ 7° =0, 
aus der hervorgeht 
tg (x — a4 = sin 4 (s) 


cos A— e448 @ 


Mittels dieser Gleichung laBt sich « als Funktion von / berechnen, woraus sich durch inverse 
Auftragung 4 als Funktion von « mit 9 als Parameter ergibt. Die auf diese Weise gewonnenen 


—————EE : 
0 30 60 90 120 490 180° 


ae 


Bild 4. Stromflu&dauer A fiir verschiedene Werte Y/R = tge, als Funktion von «. 


| 
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oe 


Bild 5. Ziindwinkel « und Léschwinkel « +4 fiir verschiedene Werte X/R = tg oe, aufgetragen tiber a. 


Kurven sind in Bild 4 wiedergegeben. Den dementsprechenden Ziindwinkel « und Lésch- 
winkel « + A zeigt Bild 5. Darin sind die. Parameterwerte zu X/R = tg 0 = 0; 0,2; 0,5; 1; 
2; 10; 0 angenommen. 

Fiir kleine Reaktanzwerte (0<<1, X < R) findet man fiir (5) die Naherung 


A=a—(a—0), (5) 


fiir groBe Reaktanzwerte >=, x> R) gilt naherungsweise 


IU 
A A ae 

Cigtd +2) 16-2 | (5b) 
= 
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3. Effektivwert des Wechselstromes 


Der Effektivwert J,, des Wechselstromes fiir « <@ hat den konstanten Wert (3). Fir 


« > ist der Effektivwert aus 
tA 


2 *R 
= = cos? @ | [sin (¥ —o) — sin (a — 9) e— @—* ee]? dx 
14 


2 


zu errechnen. Die Integration ergibt zunachst 


Wetcos 
ean 
Spaoin” (G6) tio (tee 74788) 


Mit Einfiithrung der Léschbedingung 


[A — sin A cos (2 (« —e) +A)— 4 sin (« — @) sin g (sin a — sin (a + A) e~*%8?) 


sin (a — 0) e*%@ = sin (a —o + A), 
entsprechend der Gl. (5) laBt sich der Ausdruck in der eckigen Klammer noch wesentlich 
vereinfachen. Der dritte Bestandteil in dieser Klammer wird dadurch 
= 4sinj-sing-sin (2«—o-+A) 
und der vierte Bestandteil wird 
= —sindj-tgo-sin (2(«—o)+/). 


Zieht man die drei Ausdriicke zusammen, so entsteht schlieBlich 


Ser A sin A cos (2a +0 +A) 6 
pee G pss n cos (0) 
oder, bezogen auf den Effektivwert J,,) nach (3) des ungesteuerten Stromrichters, 
fle A sin A cos (2«+ 0 +4) 
a= — (7) 
Wee M4 as cos@ 


Die Auswertung der letzteren Gleichung ergibt die in Bild 6 aufgetragenen Kurven fiir 
ld Sq Tf 4 aS U,, | R? + X2. Eine ahnliche Gleichung des Effektivwertes des Gleich- 
stromes eines entsprechend gesteuerten Gleichrichters ist bereits friither [5] abgeleitet worden. 


0 30 60 90 120 150 180° 
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Bild 6. Effektiver Strom I, bezogen auf Iq = Uy|V R? + X®*, fiir verschiedene Werte X/R = tga, als Funktion von a. 


4. Wirkleistung und Blindleistung 


Als Nachstes sollen die Wirkleistung und Blindleistung auf der Verbraucherseite berechnet 
werden. 

Auch diese Rechnungen enthalten im Endergebnis e-Funktionen, die sich mit Anwendung 
der Léschbedingung (5) auf Kreisfunktionen zuriickftihren lassen. So findet man fiir die 
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Wirkleistung 


ata 
2 
Py = aoe COS @ | [sin (v — 9) — sin (a — 9) e@—8e] sin x dx 


aw 


2 
=e [A cos@ —sinAcos (20 +0+A)], 


welcher Ausdruck die Beziehung Pjy = R Ji, mit J,, nach Gl. (6) in sich schlieBt, und fiir 
die Blindleistung 
5 ata 


U : J 
a mas Cos @ | [sin (¥ — e) — sin (a — 9) e*—*)“£2] cos x dx 


a 


= T2901, sin g —sin Asin (20 -+e+)]. 


cose = Oe ao 0 Sd 


ist, wenn J,,9 den effektiven Wechselstrom fiir « = 0 bis « = @ und Ps» die Scheinleistung 
des Verbrauchers fiir « = 0 bis « =@ (d.h. die ungesteuerte Scheinleistung) bedeutet, so 
1aBt sich schreiben 

) 


P : 

MEAL eons ap eda ee pe A) ; 

Pso I 5 (8) 
Pp A si 


—_- —— sin 9 — - 
Pso wt e 


24 sin (20 -+@ +A). 


Diese GréBen findet man in Bild 7 und Bild 8 kurvenmaSig aufgetragen. 


= —— # 1 ! 1 | 
0 30 60 90 120 150 180 0 30 60 99 120 150 
C—- ——— 
Bild 7. Wirkleistung Py, bezogen auf die Scheinleistung Ps 9, Bild 8. Blindleistung Py, bezogen auf die Scheinleistung Ps 0 
fiir verschiedene Werte X/R = tgo, als Funktion von a. fiir verschiedene Werte X/R = tgg, als Funktion von «. 


In allen Kurvenblattern ist eine gestrichelt gezeichnete Grenzlinie eingetragen, die den 
Beginn der Aussteuerung bei a =e anzeigt. Fiir « <@ sind die Funktionswerte konstant 
und haben den fiir den ungesteuerten Stromrichter geltenden Wert. 


5. Anoden-Sperrbeanspruchung 


Die in jedem Léschaugenblick der Ventile auftretende Anoden-Sperrbeanspruchung ist 
nach der von Wasserrab gegebenen Definition bestimmt durch [6] 


‘dt 
B=AUsy+(S*) (0) 


180° 
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worin AUg, die im Léschaugenblick erscheinende Sprungspannung ist. Da allgemein 


X diy ; : 
—— + Kt, = 72 0, sin oF 
w dt q Yo / w oil a 
und im Léschaugenblick 7,, = 0 ist, so ist AUs, = V2 U, sin (a +4) =X (FI, -. 4 und 
damit F 
4)? U%,+ ; 
B= (4 Usp) oz — ‘ pa sin? (ov + A) . (10) 


xX x 


Die Gleichung B = (AU<5,)? w/X = (AUs,)*/L ist iibrigens eine ganz allgemein geltende Be- 
zichung, die fiir alle Stromrichter mit einer Induktivitat in den Ventilleitungen zutrifft. 


20 


45 


1 
0 30 60 90 120 150 180° 
Bild 9. Anoden-Sperrbeanspruchung B, bezogen auf w Ps 0 fiir verschiedene Werte X/K = tg oe, als Funktion von a. 


Da mit tgo = X/R 
Uw = Up Iwo Ee Pso 
ary sin 0 sin 0 
ist, wird die Anoden-Sperrbeanspruchung, bezogen auf die ungesteuerte Scheinleistung P.y, 


Boy 2 sink ie A) (11) 


wo Psy sin 0 
Diese Beziehung ist in Bild 9 aufgetragen. 
Mit Verwendung der ersten Gleichung von (8) ist es ohne weiteres méglich, die Anoden- 
Sperrbeanspruchung beispielsweise auf die Wirkleistung zu beziehen. Es wird dann 


ely, Pee (ee (12) 
ore open sin A ; 
sin 0 (= cos 0 — - cos (24 + A+ 0) 
Lh [ 


Damit sind alle GréBen, die fiir die Auslegung des Stromrichters von Bedeutung sind, fiir den 
praktischen Gebrauch aus den gegebenen Kurven bestimmbar. 

Bei Verlust des Sperrvermégens der Ventile tritt keine Riickziindung im iiblichen Sinne 
mit nachfolgendem TransformatorkurzschluB ein, sondern nur ein Vorspringen der Ziindung 
und damit ein Anstieg des Wechselstromes /,, auf den ungesteuerten Wert Jig - 
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Die wechselstromseitigen Spannungsharmonischen eines 
Stromrichters bei idealer Sperresonanz des Netzes 


Von 
DIETER RuMPEL, Erlangen 


Mit 12 Textabbildungen 


(Eingegangen am 30. Juli 1960) 


Ubersicht: Es wird berechnet, welche Héhe die Spannungsharmonische auf der Wechselstromseite eines 
gesteuerten Stromrichters erreicht, wenn die zugehdérige Stromharmonische durch eine Netzresonanz ideal 
gesperrt ist. 

Zur Vereinfachung der Rechnung wird angenommen, da8 das Netz bei der jeweils betrachteten harmoni- 
schen Frequenz eine unendlich hohe Sperrimpedanz aufweist, wahrend es sich fiir alle anderen Frequenzen 
wie eine Induktivitat verhalt. Fiir den Stromrichter wurde vollkommene Drosselglattung des Gleichstroms 
und Verlustfreiheit der Ventile vorausgesetzt. 


Ubersicht der benutzten Formelzeichen 


eee Maximal-, Scheitelwert eines Stromes 

ib, Induktivitat 

Oh, u Effektiv-, Scheitelwert einer Spannung 

Dy; Dy, FourtERkKoeffizienten der Ordnungszahl p 

a Augenblickswert eines Stromes 

iy Harmonischer Stromanteil der Ordnungszahl p 

m Zahlindex der Ventile m = 0... (p— 1) 

nN beliebige ganze Zahl 

p Pulszahl des Stromrichters 

u Uberlappungswinkel 

% Zeit, gemessen im Winkelmaf der Grundschwingung 

Xr PRr YR Hilfsfunktionen 

€ Zindwinkel (im allgemeinen mit « bezeichnet) 

lL Sperrfrequenz des Netzes bezogen auf Grundfrequenz 
v Ordnungszahl der Harmonischen 

P Schwenkwinkel des Bezugsventils 

Index g Gleichstrom, kommutierter Gleichstrom 

Index k kommutierender Grundschwingungs- KurzschluBstrom 
Index w kommutierender Strom von Sperrfrequenz 


Index Rk, S, T Kennzeichen der Drehstromstrange 


Die GréB8en in den Ventilkreisen sind einfach, die Gr6éBen in den Netzkreisen zweifach gestrichen. 


1. Einleitung 


Bei der Projektierung groBer Stromrichteranlagen taucht haufig die Frage auf, welche 
Riickwirkungen die im Betrieb auftretenden Stromharmonischen auf den Netzspannungs- 
verlauf ausiiben, insbesondere, wenn das Netz bei der Frequenz einer der Harmonischen eine 
Sperresonanz aufweist. Es sind deshalb wiederholt Arbeiten veréffentlicht worden, die rechne- 
rische oder graphische Verfahren zur Ermittlung der Wechselspannungsharmonischen ent- 
halten. 

So wird in zwei amerikanischen [1, 2] und einer russischen [6] Verdffentlichung vor- 
geschlagen, zur Ermittlung der Harmonischen den Stromrichter fiir jede auftretende Ober- 
frequenz getrennt durch einen aktiven Zweipol zu ersetzen. Die Innenimpedanz dieses Zwei- 
pols ist ohmisch und seine Leerlaufspannung (bezogen auf die grundfrequente Spannung) 
entspricht gréBenordnungsmaBig dem auf Grundschwingungsstrom bezogenen Strom der betr. 
Harmonischen bei idealer Kommutierung. Wenn das Netz bei der betrachteten Oberfrequenz 
eine Sperresonanz aufweist, so wiirde dieser Zweipol cine Klemmenspannung von max. der 
GréBe seiner Leerlaufspannung liefern. Das bedeutet, da8 beispielsweise im ungiinstigsten 
Falle eine 11. Spannungsharmonische von 100/11 ~ 9% im Netz zu erwarten ware. 
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Andererseits wird von LEUKERT und KUBLER [3] vorgeschlagen, die harmonischen Strom- 
anteile des Stromrichters als eingepragte Stréme aufzufassen, so daB je nach Eingangsimpedanz 
des Netzes bei den betreffenden Oberfrequenzen Spannungsharmonische beliebiger Héhe ent- 
stehen kénnten. 

Auch AyMANNS [4] nimmt an, daB die ideale Sperresonanz des Netzes bei einer Ober- 
frequenz durch das Auftreten einer unendlich hohen Spannungsharmonischen gekennzeichnet 
ist; obwohl dieses Verhalten nicht notwendig aus 
seinem Rechenansatz hervorgeht. 

In einer Arbeit von KUBLER [5] wird die Hohe 
der netzseitigen Spannungsharmonischen unter- 
sucht, wenn man einem Stromrichter einen rein 
sinusférmigen Netzstrom aufzwingt, d.h. alle 
oberfrequenten Stréme gleichzeitig sperrt. Aus 
diesem Ansatz folgt notwendig, daB jeweils zwei 
Ventile gleichzeitig Strom fithren und der Gleich- 
strom einen welligen Verlauf aufweist, der vom 
Netzstrom fest vorgegeben ist; der Einsatz von 
Glattungsdrosseln fiihrt lediglich zur (beliebigen) 
Erhéhung der netzseitigen Spannungsharmoni- 
schen. Auch dieses Ergebnis — welches naherungs- 


\z\ 


Bild 1. Frequenzgang des Betrages der vor dem Stromrichter 


ee eee weise das Betriebsverhalten eines Stromrichters 
bei Sperrung seines starksten Stromharmonischen- 
paares wiedergibt — weist darauf hin, daB die Spannungsharmonischen bei Netzresonanz 


u. U. betrachtliche Werte annehmen kénnen. 

Die vorliegende Arbeit stellt einen Ausschnitt aus der Dissertation des Verfassers dar. 
Es wird versucht, auf Grund eines weitergehenden Ansatzes die Spannungsharmonischen bei 
Sperresonanz zu ermitteln, wobei insbesondere die Steuerung des Stromrichters und der Ein- 
flu8 der Spannungsharmonischen auf den Kommutierungsvorgang beriicksichtigt werden. 
Natiirlich muBten die auBerst komplexen Zusammenhange auch hier z. T. stark idealisiert 
werden, um den Rechenaufwand in ertraglichen Grenzen zu halten. Die Idealisierungen 
betreffen vor allem den Frequenzgang des Wechselstromnetzes, wo vorausgesetzt wurde, daB 
das Netz einschlieBlich Stromrichtertransformator bei der jeweils betrachteten harmonischen 
Frequenz eine unendlich hohe Sperrimpendanz aufweist, wahrend es sich fiir alle anderen 
harmonischen Frequenzen wie eine Induktivitaét verhalt (Bild 1). Fiir den Stromrichter 
wurde vollkommene Drosselglattung des Gleichstroms und Verlustfreiheit der Ventile voraus- 
gesetzt. 


2. Betrachtete Stromrichterschaltung 


Bild 2 zeigt schematisch einen Stromrichter in p-pulsiger Sternpunktschaltung, der iiber 
die Induktivitaten L” an ein starres Drehstromnetz (U,, Us, U) angeschlossen ist. Die prak- 
tische Ausfiihrung des Strom- 
richtertransformators hat auf 
das Ergebnis der folgenden Uber- 
legungen keinen Einflu8; der 
Einfachheit halber sei angenom- 
men, daB er in Form eines 
Drehtransformators mit dem 
Ubersetzungsverhiltnis 1:1 aus- 
gefiihtt ist, bei dem die ventil- 
seitigen Wicklungen symme- 
trisch tiber den Umfang verteilt 
sind. Die Induktivitat L’ enthalte alle im Kommutierungskreis wirksamen Induktivitaten. 
Das Filter, welches in Verbindung mit der Induktivitat L” den idealisierten Frequenz- 


Bild 2. p-pulsiger Stromrichter in Sternpunktschaltung. 
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gang nach Bild 1 realisiert, 1aBt sich als dampfungsloser Sperrkreis auffassen, dessen Induk- 
tivitét gegen Null und dessen Kapazitat gegen Unendlich gegangen ist derart, daB das 
Produkt beider GréBen konstant und damit die Sperrfrequenz erhalten blieb. (Man kann das 
Filter aber ebensogut als eine eingepragte konstante Spannung Uy von der Frequenz der zu 
sperrenden Harmonischen ansehen. Gro8e und Phasenlage der Spannung sollen dann so be- 
messen sein, daB der von dieser Spannung 
getriebeneStrom gerade alle harmonischen 
Stromanteile der Sperrfrequenz im Filter- 
zweig zu Null erganzt). Da — wie in der 
Rechnung angenommen — im Gleich- 
stromkreis eine unendlich hohe Glattungs- 
induktivitat liegt, kann die Spannung am 
Filter nur wahrend der Uberlappungszei- 
ten einen Strom treiben, der sich zum kom- 
mutierenden Grundschwingungsstrom ad- 
diert und somit auch die Dauer der Kom- | z 

| 


mutierung beeinfluBt. ae 


Der kommutierende Strom von Sperr- . 
frequenz wird — ebenso wie der kom- ay 
mutierende Grundschwingungsstrom — ; 
von der Induktivitaét ZL” begrenzt und 
weist ein vom Ziindzeitpunkt abhangiges 
Gleichstromglied auf. Auf Grund der 
Kommutierungsbedingungen des Strom- 
richters ist es méglich, den Verlauf dieses 
Stromes bei idealer Sperrung zu berech- | 
nen. Damit ist gleichzeitig die Héhe der | 
treibenden Spannung U; festgelegt, welche | 
im folgenden kurz als ,,Sperrkreisspan- | 


nung‘‘ bezeichnet werden soll. yp 
vr \ 
Vv Vv 


hal 
3. Die Stromharmonischen | | aa 
in den Ventilkreisen a | 
Bie 


Ue 


Zunachst wird der Stromverlauf auf 
der Ventilseite und hier wiederum in n | \ ton 
einem ,,Bezugsventil (m= 0) betrach- 
tet. Bild 3 zeigt oben den Verlauf des Bild 3. Stromverlauf im Bezugsventil (#,.= 0) bei Sperrung 
Ventilstromes 7’ bei Sperrung einer Har- elas Heroes 
monischen, wobei als Abszisse die Zeit 


und als Koordinatenursprung das Ende der Kommutierung gewahlt wurde. Dieser Strom- 
verlauf setzt sich zusammen aus: 


dem kommutierten Gleichstrom 7, 

dem kommutierenden Grundschwingungs-KurzschluBstrom 7; 

dem kommutierenden Strom 7, von der Frequenz der gesperrten Harmonischen v = y, 
der von der Sperrkreisspannung U; getrieben wird. 


ZweckmaBigerweise sind die kommutierenden Teilstréme 7, und 7, am Ende (I) und Anfang 
(II) der Brenndauer in Bild 3 getrennt aufgetragen. 


Der Ventilstrom 7’ enthalt einen harmonischen Anteil 7; von Sperifrequenz (vy = y), der 
sich mit seinen FourtERkoeffizienten a, und 0; schreiben 1aBt: 


i, = a,cosyx + b,sinvx. (1) 
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Die Fourterkoeffizienten ai, 6; des gesamten Ventilstromes setzen sich additiv aus den FOURIER- 
koeffizienten der Teilstréme 7, 7, 7, fiir die betr. Harmonische zusammen; man kann fiir (1) 
schreiben 
ty = (az, + iy + ay,) cosy x + (b,, + by, + b,,) Sin Vx. (2) 
Die Fourterkoeffizienten der Teilstréme findet man wie folgt: 
Kommutierter Gleichstrom: Aus dem Verlauf des Stromes 


A 2 mu 
~Z2= J, OS ee at tee (3) 
liefert die FourtERanalyse die Koeffizienten: 
ts Ft gin ?® =2sin™. J cog ™ , 
YI Pp p VU (4) 
iy Je 20 Se NS Peat | 
gy = ——*£]1 — cos 2sin—.- —£sin—. 
2 VI P pb var 
Kommutierender Grundschwingungsstrom: Aus dem Verlauf des Stromes 
i,, = Ji, [cos (v + u + 2) —cose] FSU TIC; (5) 
liefert die FouRIERanalyse die Koeffizienten: 
ay,—Tha, tir = 4B, (6) 
14 4 


a, und f, sind zwei der spater behandelten Hilfsfunktionen. Der Teilstrom 7,,;; geht aus 7, 


durch Verschiebung um“ und Wechsel des Vorzeichens hervor; fiir seine FourrERkoeffi- 
zienten gilt daher 


a 
/ 


a 25 : 2% 
ee a, cos +B siny =), 
4 : ; (7) 
csp 
4 PAGE a 270 
bj = —J8)g cosy +a sine 
kIIv k k 
‘3 p 2 


Damit kann man zusammenfassen: 
ey 
Jk 


| 


‘ 
aky 


/ , : M4 
Agty + Ary = 2 Senses , 


: M4 4 
a, Sin » — — B, cos »— 
P P 


| 


Dey = Dery + 0; ,=2sinyS . sin y= a,cosr 2]. 
k kI RIT p e B, res k p 


Kommutierender Strom von Sperrfrequenz: Von diesem kommutierenden Strom 
ist zunachst nur bekannt, daB seine Frequenz der Ordnungszahl der gesperrten Harmonischen 


/ entspricht und daB er ein Gleichstromglied aufweist, welches seinen Anfangswert bei x = — u 
zu Null erganzt. Der allgemeine Ansatz hierfiir lautet: 
tor =C' + D'cosux+E' sinus %=—U...0, (9) 
mit dem Gleichstromglied C’ = — D' coswu-+ E'sinwu. 
Die FourtERanalyse liefert bei diesem Strom fiir y = w die Koeffizienten: 
Got» = = _[v u—cos» 4+ sin vu] + ——sin? » a, | 
2V7 2VIN 


(10) 


(3 e 


[1 + cos? yu —2cosvu] + 
21 2V7 


1 


Corrs = 


[vy uw — (2— cos v4) sin u] , | 


t,71 geht aus 7, durch Verschiebung um = und Wechsel des Vorzeichens hervor [7]; fiir 


den harmonischen Stromanteil von Sperrfrequenz 7,7; lauten daher die Fourterkoeffizienten 


‘ 2 2% , . 27 
Cant ty = a Gare cas Rear lores > 


(11) 


/ 


b = , 20 , roe 
ON A OS ae bath COS ee wT an eRe 
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Als Summe der Fourierkoeffizienten (10) und (11) erhalt man schlieBlich 


, , ‘ It 
Qe» = Cole + ay lly = ae os 


, . 4 jy MLA 
Ay tySin ¥ —— by 7, COosv | | 
Pp p (12) 


, , + MA ° A , MA 
Oe 0574+ lat is = sine 2 [0 1,sine = + digs, 0089 7 | | 


Mit diesen Fourterkoeffizienten der Teilstréme ist der Verlauf der Stromharmonischen von 
Sperrfrequenz im Bezugsventil (m = 0) festgelegt. 


Das m-te Ventil arbeitet — um x = m = versetzt — spater. Es folgt also: 
ism = a, cos (x—m 22) + 6, sin »(x— m2”) (13) 
4. Hilfsfunktionen 
Die Fourierkoeffizienten und der Verlauf des kommutierenden Kurzschlu8stromes 


hangen proportional von dessen Scheitelwerte vis und nach verhaltnismaBig komplizierten 
transzendenten Funktionen vom Steuerwinkel ¢ und vom Uberlappungswinkel u ab. Fiir die 
praktische Rechnung wurden daher folgende Hilfsfunktionen eingefiihrt, die in [7] tabelliert sind. 


1. Die bezogenen Fourterkoeffizienten «, und f, (vergl. (6)) 


hee sin (u + ¢) sin (vu + 8) sin (v u— €) 
i 4 200 Lae 2964), 7 (14) 
v cos (u + &) cos (v u + €) cos (vy u— €) Cos € 
ica ye —— 4 15 2v(v—1) gente? v | 


2. Der auf ft bezogene Momentanwert, den der kommutierende KurzschluBstrom am 
Ende der Kommutierung (x = 0) aufweist 


(ie) 


v= ae cos (E+ 2): — COs &. (15) 
Ji ; 
Nahert sich — wie es bei gréBerem Steuerwinkel und kleinem Uberlappungswinkel der 


Fall ist —— der Verlauf des kommutierenden KurzschluBstromes einer Geraden, so kénnen 
fiir (14) auch folgende Naherungsformeln benutzt werden: “ 
yR 


a, & mad [1—cosvu], | 
vu 
Vk ; 
w— —* [yu—sinr 4]. 
Bi ee Age, ‘ | 


sit rit 
ly yg 


5. Summation der Ventilstromharmonischen 
zur Netzstromharmonischen 


Die Stromharmonischen in den Ventilkreisen sind — im Ge 
gensatz zu den eigentlichen Ventilstrémen, aus denen sie er 
mittelt wurden — Dauerstréme und kénnen daher wie normal® 
einfrequente Wechselstr6éme behandelt werden. Bei der ge- 
wahlten Ausfithrungsform des Stromrichtertransformators ist m 
das Ubertragungsverhalten zudem besonders anschaulich, ohne ae 
da8 die Allgemeingiiltigkeit des Ergebnisses beeintrachtigt Derchihtnesschees alee Witla 
wiirde. Bild 4 zeigt das raumliche Schema der Durchflutungs 
achsen der Transformatorwicklungen. Bei der Ubertragung eines Ventilstromes #,,, in den 
Primarkreis miissen zwei Bedingungen erfiillt werden: 


1. Die Summe der Primarstréme ist Null 


° St ss 
Letts +t7r=0. (17) 
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2. Die geometrische Summe der Durchflutungen ist Null 


4, R COS (p +m “2 + 47, cos (r yy ee 2) + 4) 7 cos (p VES gence =| Sater 
p p 3 P 3: (18) 
tp sin (p +m 3) go sin (p +m a ae ‘2 Sy Re sh (r +m = S =O. 


Aus diesen Bedingungen erhalt man fiir die aus dem m-ten Ventil in die Primarstrange 


tibertragenen Stréme 


“7 RSC 27 

VR 3m cos |p + ™ A : 

wt Th 47 2 IU 

1s 5 — dy, COS (r +7" 4+ m — ; (19) 
3 3 P 

wt Dey ZH ZI 

ae 7 i,m 008 (p+ Sri Meenas 


Aus (19) folgt fiir den aus 
Strom der v-ten Harmonischen: 
(pen 


allen Ventilen beispielsweise in den Leiter R tibertragenen 


al 2 ” 2H 
Deri te tym COS m —— 
: 3 m=0 G # | 
und mit (13) 
p-1 
“ Zz , 27 fs 27 2 
a a, COSY |X% — m — b, sin v x — m=) cos ( m 22) 20 
ser eer ee mig 
Die Auswertung dieser Summe fiihrt auf 
ifn = © {a, cos (vx tq) +0,sin (vx + 9)}, (21) 
fiir 


n ist hierbei eine beliebige ganze positive Zahl; die Vorzeichen in (21) und (22) miissen iiberein- 


stimmend gewahlt werden. 


Fir alle Ordnungszahlen auBer (22) ist die Summe (20) gleich Null; Harmonische dieser 
Ordnungszahlen treten bei der betrachteten Stromrichterschaltung nicht im Drehstrom-Netz 


auf. 


Entsprechend findet man fiir die Harmonischen in den beiden anderen Primarstrangen: 


oS = P lal 005 v% +(9 + *2) +- b, sin 
ifn =? la, cos vx + (p+) oh Sin 


v3) 


aa) 


6. Ableitung der Bestimmungs¢gleichungen 


Fiir die Stréme nach (21) und (23) mu8 entsprechend der Schaltung von Bild 2 bei » = w 


gelten: 


te =ts=t47=0. 
Damit erhalt man aus (21) und (23) die Bedingung 
he = 0,0), 
oder mit (2) 
gy + Gy + by, = 0; 
bey + Ory + On, = 0. 


(24) 


(25) 


(26) 
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Durch Einsetzen der FourtERkKoeffizienten fiir die Teilstréme (4), (8), (12) ergibt sich: 


IE : MA 2 Jt 2 ery 

sin »y — (72 ay + 2 dp1,)-sin » = — ay te Je | 2 b,,1,\Cos vy —| =0, 
p % p a v1 p 

siny™ ee ele sin» ~ 2Jh +24, Zs sk 
1 | Ls aE : Dee hee ty)? cosy —| =0. 
p a k Poe wily p = k Sha b 


Die Gleichungen (27) sind wegen (22) fiir die auftretenden Harmonischen nur erfiillbar 
durch 


2d , 
a Ce 2 Oy py 05 | 
~ (28) 
/ 
PD Vlys 2] \ / 
: 2051, = 0. 
ze By aoe | wily 
Setzt man hierin die Gleichungen (10) ein, so erhalt man schlieBlich: 
2 Ji D’ ; Ey 
: X, + —[vyu—cosvusiny u] + —sintryu=o, 
4 Vv VI ( ) 
a , 29 
Zack D in ; 
2Jag Se — -— [1 + cos?» 4 — 2 cosy u] + — [vy u— (2—cosvu) sinvyu] =o. 
5 vn va v7 


Als weitere Beziehung findet man an Hand von Bild 3: Die kommutierenden Teilstréme 
i, und 1; miissen am Ende der Kommutierung (« = 0) gerade J, zu Null ergianzen, d. h. 
Vode Dis Gea= FeO. (30) 

Bei einer bestimmten Stromrichteranordnung sind in der Regel der Gleichstrom J,, der 
Scheitelwert des kommutierenden Grundwellenstromes ii , und der Steuerwinkel ¢ gegeben; 
gesucht sind die simultanen Werte von D’, E’ und der Uberlappungswinkel w bei idealer 
Sperrung einer Harmonischen. Grundsatzlich sind die gesuchten Werte durch die 3 Glei- 
chungen (29) und (30) festgelegt; die praktische Berechnung st68t jedoch auf erhebliche 
Schwierigkeiten, da sich der Uberlappungswinkel w nicht explizieren 1aBt. 

Man geht deshalb zweckmaBig so vor, da8 man # und ¢ als bekannt voraussetzt und die 
zugehérigen Werte von eure » Def f. und E’/ ie ermittelt. Fiir die tabellarische Berechnung 
der Beziehungen ist diese Vertauschung von gesuchten und bekannten GréBen ohne Belang. 
Aus (29) und (30) ergeben sich so die folgenden Bestimmungsgleichungen 


Jeg _ 2a, (1— cosy u) + 2” B, siny hemes, (vy u— sin v u) 


We yu-+sinvu 


9 & - 
—2va,[yu—sinvu- (2—cosv u)] + ( v Br—2 ee sint 
k 


(vy w—sin vu)? 


2 Jeg i 
— 2” a, [1 + cos* » wu— 2c0s vy u]— {2 Bp — 2 i, {vy w— cos? u sin » u] 
k 


- j i , 


Fe (vy u— sin v u)* 


woraus man unter Verwendung der Hilfsfunktionen verhaltnismaBig leicht die Quotienten 

, Be 

= und x (32) 

Jk Jk 

remitteln kann. Diese beiden Werte legen den gesuchten Verlauf des kommutierenden Stromes 

von Sperrfrequenz fest, der die Stromharmonische von Sperrfrequenz im Filterzweig zum Ver- 

schwinden bringt. 

Es besteht nun noch die Aufgabe, aus diesen Werten auf die Sperrkreisspannung U, zu 

schlieBen, die den Strom 1, treibt. 
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Mit den ermittelten Quotienten (32) erhalt man 
RSS 


Dp’ \2 2 
Vey eh. 
Jk ie Jk 


worin Te der Scheitelwert des Wechselstromes in 7, ist. Die beiden Stréme 7, und 4, flieBen 
zu gleichen Zeiten in den gleichen Strombahnen und werden beide durch L” begrenzt. Es 
muB also gelten: 


(33) 


as 
Ge wio lt? ee: 
ep ee a" — a” (34) 
Jk tes y Uy 

w ie 


_Fiir die auf Grundschwingungsspannung bezogene Spannungsharmonische bei idealer 


Sperrung findet man damit: 
ut WeNe cm Wee 2 
Vi Gy 
Uy Jk Jk 


Das Dach iiber den Spannungssymbolen kann hierbei auch fortgelassen werden, da der 
Quotient aus den Effektivwerten der gleiche ist. 


> 


= 


(35) 


7. Die ermittelten Sperrkreisspannungen 


Mit Hilfe der Gl. (31) und (35) wurden die Sperrkreisspannungen fiir die ungeraden Ord- 


nungszahlen zwischen 3 u. 13 in Abhangigkeit von den Parametern ia 'J;, und € berechnet. 
Sie sind in den Kurvenblattern Bild 5—10 aufgetragen. Die Bilder enthalten weiterhin den 
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Bild 5. Sperrkreisspannung fiir die 3. Harmonische. Bild 6, Sperrkreisspannung fiir die 5, Harmonische. 


i meee 402 Oe eh = *-60° 80° 
9 
; ja 
3010 
6 50% 
5 70 
» 80 
4 a 90 
5 100 
R) gee 
110 
2 
420 
BS 
=4 
40 
: 425 
é 121% 
5t—-S 
4 105 
46 
Z 
5 % 
2 2 Ks 
10. —|—____ 
9 r) 
8 2 
& Be 
5 of 
4, Vs | 
3 
2 
40° at 
Bild 7. Sperrkreisspannung fiir die 7, Harmonische 
-5° 0° “0pm eDveeres0e a 4o° 80° 60° e 
4 80 
g 
é | 
: 30% 20% 
5 S LES 
Bel SE 10 
f Clo 90 
iy 4 
2 00 
gs 10 
f 0 
—4 
10, 
7 » 
H | 
5 
X [ 
3 T 
2. 
: 2 
css 
2 S 
0 ; = = == 
Je 
6 = 
5 s 
4 
3 
2 T Hs 
° 


Bild 9. Sperrkreisspannung fiir die 11. Harmonische. 
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Bild 10, Sperrkreisspannung fiir die 13. Harmonische. 
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bei idealer Sperrung der betr. Harmonischen zu beiden Parametern gehérigen Uberlappungs- 
winkel wu. 

Die errechneten bezogenen Sperrkreisspannungswerte gelten fiir alle Stromrichterschal- 
tungen, unabhangig von der Pulszahl und der primaren Phasenzahl [7]. Bei Saugdrossel- und 
Schwenkschaltungen gelten sie allerdings nur solange, wie die einzelnen Kommutierungs- 
gruppen gleich belastet sind und sich nicht gegenseitig in ihrer Kommutierung beeinflussen. 

Die effektive Pulszahl bestimmt nach (22) nur, welche Ordnungszahlen iiberhaupt auf- 
treten, nicht aber die Héhe der Sperrkreisspannung fiir die auftretenden Harmonischen. 


Aus der Betrachtung der Kurvenblatter ergeben sich 5 wichtige Folgerungen: 


1. Die Sperrspannung einer Harmonischen ist weitgehend abhangig vom Verhdltnis Jide 
d.h. vom Verhaltnis der Stromrichterleistung zur Grundschwingungs-KurzschluBleistung im 
Kommutierungskreis. 


2. Die Sperrkreisspannung steigt mit wachsender Ziindverzégerung stark an. 


3. Der Uberlappungswinkel ist bei Sperrung einer Harmonischen zum Teil wesentlich 
gréBer als bei induktivem Netz. 


4. Bei kleinem Ziindwinkel sinkt die Sperrkreisspannung mit zunehmender Ordnungszahl 
der Harmonischen ab, wahrend sie bei groBem Ziindwinkel ihren Wert etwa beibehalt. Z. B. 


fir Je =O nll 


Jk 
« | user | vgjor | vzsor | usr | enter | vyluy 
0° 73% 50% 37% 27% 19% 16% 
30° | 106% 96% 88% 81% 73% 69% 
60° 124% 122% Oe 118% | 115% 112%, 


5. Die Sperrkreisspannung weist fiir 
vu Fw 360° , 


sehr geringe Werte auf; auch der Uberlappungswinkel ist in diesem Gebiet bei Sperrung nur 
unwesentlich gegeniiber dem induktiven Fall vergr6Bert. 


8. Der Kommutierungsbeginn bei Sperrung einer Harmonischen 


Bei der Berechnung der Sperrkreisspannung wurde angenommen, daB der Kommutierungs- 
beginn (x = — ) willkiirlich durch die Steuerorgane des Stromrichters festgelegt ist. Bei 
Sperressonanz ist jedoch der Spannungsverlauf an den Ventilen durch die iiberlagerte Span- 
nungsharmonische erheblich verzerrt; es muB daher noch der Gesichtspunkt beachtet werden, 
daB die Kommutierung — unabhangig von der Freigabe des Ventils durch die Steuerung — 
nur unter der Bedingung einsetzen kann, daB die Spannung an dem neu zu ziindenden Ventil 
in DurchlaBrichtung hegt. Die aquivalente Bedingung fiir den kommutierenden Strom lautet 
mit den hier eingefiihrten Gré8en: Der kommutierende Gesamtstrom (7,; + ¢,;) kann nicht 
mit positiver Tangente beginnen; mathematisch gefaBt 


dik diy 
(eae! GE | eee So, (36) 
oder mit (5) und (9) 
—J,sine +yD'sinvu+tvE' cosvu<o. (37) 


Bei kleinen Werten von ¢ ist diese Beziehung nicht immer gegeben; in solchen Fallen folgt 
der Kommutierungseinsatz nicht mehr der Steuerung, sondern bleibt bei dem gréBeren Winkel 
e* liegen, fiir den gilt: 


LTA / : ¢ 
— J, sine* + 7 Dis siny u + » Eis) cosyu =o. (38) 
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Ist der Steuerimpuls zu kurz, so daB fiir e* das Gitter bereits wieder sperrt, so bleibt die 
Kommutierung tiberhaupt aus. Beim ungesteuerten Gleichrichter beginnt die Kommutierung 
immer bei € = é*. 

Bild 11 zeigt den ungefahren Verlauf von e* bei idealer Sperrung einer der untersuchten 
Harmonischen — soweit sich dieser Verlauf aus den Werten ermitteln lieB, die bei der B>- 

ee rechnung der Sperrkreisspannung anfielen. Bei der Sper- 
4 40°_-30" _-20" _-10" __o?__9°_2°- ss rung der Harmonischen niederer Ordnungszahl liegt «* 
durchweg im negativen Gebiet. Der Kommutierungseinsatz 
folgt hier der Steuerung sogar bis zu negativen Ziind- 
winkeln, beim ungesteuerten Gleichrichter setzt die Kom- 
mutierung infolge der Sperresonanz verfriiht ein. 


9. Experimentelle Uberpriifung der Rechnung 


Zur Uberpriifung der Rechnung wurde zwischen einen 
Zweipulsgleichrichter kleiner Leistung und das Netz eine 
variable Induktivitat und ein auf die 5. Harmonische ab- 
gestimmter Sperrkreis geschaltet. Um eine méglichst gute 
Annaherung des in Bild 1 gezeigten idealisierten Frequenz- 
gangs zu erhalten, wurde die Sperrkreisinduktivitat sehr 
klein gegeniiber der in Reihe geschalteten Induktivitat ge- 


— Ge gut belegt gusgezogene Kurven =berechnete Werte 


=Mebwerte 
——— Kurven schlecht bele 
FNS FAG 0 015 02 0,25 03 
no 
Lk 
Bild 11. Verlauf von e* bei vollkommener Bild 12. Verlauf der Sperrkreisspannung der 5. Harmonischen und des Uberlappungs- 
Sperrung einer Harmonischen. winkels bei e = 0 als Funktion von IglF;,. 
& 


wahlt; durch eine elektronische Schaltung wurde der Sperrkreis weitgehend entdéampft. Durch 


Variation der Reiheninduktivitat konnten verschiedene Werte des Quotienten i, el J 4 eingestellt 
werden. 


Bild 12 zeigt die MeBwerte einer Versuchsreihe mit e = o. Aufgetragen sind der Verlauf 
der Sperrkreisspannung und des Uberlappungswinkels tiber J;,/J;,; die ausgezogenen Kurven 
stellen die berechneten Werte dar. 


10. Zusammenfassung 


Es wurde gezeigt, welche Hohe die Spannungsharmonische auf der Wechselstromseite 
einer Stromrichterschaltung erreicht, wenn die zugehérige Stromharmonische durch eine 
Netzresonanz vollig gesperrt wird. Zur Vereinfachung der Rechnung wurde ein idealisierter 
Frequenzgang des Netzes, ideale Drosselglattung des Gleichstroms und Verlustfreiheit der 
Ventile angenommen. Die Steuerung des Stromrichters sowie der Einflu8 der Spannungs- 
harmonischen auf den Kommutierungsverlauf wurde beriicksichtigt. Es zeigte sich, daB der 
auftretende Wert der Spannungsharmonischen stark vom Verhaltnis der Stromrichterleistung 
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zur Grundschwingungs-KurzschluBleistung im Kommutierungskreis, vom Ziindwinkel und von 
der Ordnungszahl der gesperrten Harmonischen abhangt. Die Spannungsharmonischen 
gehen auch bei idealer Sperrung der zugehérigen Stromharmonischen und idealer Drossel- 
glattung des Gleichstromes nicht gegen Unendlich, sie kénnen aber immerhin bei groBem 
Ziindwinkel selbst bei hdheren Ordnungszahlen noch Werte annehmen, die die Grund- 
schwingungsspannung iibersteigen. 


Literatur 


{1] Bryz, Kent: Effects of Rectifiers on System Wave Shape. AIEE Transact. (El. Eng.) 1934, S. 54—61. 


[2] Joint Subcomm. of Dev. and Res.: Wave Shape Charakteristics of Power Circuits Supplying Rectifiers. 
Engn. Rep. 22 of Joint Subcomm. of Dev. and Res. 1932, S. 183—204. 


[3] LeuKERT u. KUBLER: Oberwellenbelastung von Drehstromnetzen. E. u. M. 1936, S. 37—44, 52—55. 


{4] Aymanns: Der EinfluB der Netzkapazitaten auf die Stromoberwellen in Gleichrichteranlagen. ETZ 1937 
S. 499—503, 535—538. 


[5] Ktsrer: Stromrichterbelastung von Generatoren und Drehstromnetzen in vektorieller Darstellung. 
W.V.S. 18 (1936) S. 50—72. 


[6] Orx1na: Higher Harmonics in a Power System Supplying Mercury Rectifiers. Direct Current (1955) 
S. 115—121 — [Elektritschestwo. Nr. 2 (1955)]. 


[7] RumpeEL: Dissertation, Technische Universitat Berlin, 1959. 


Dr.-Ing. DrETER RuMPEL, in Fa. Siemens-Schuckertwerke, Technische Stammabteilung, Erlangen 


<LVI, : ; 
Hen aeks B, Berke’: Lebensdauer eines Wolframheizfadens 59 


Die Lebensdauer eines Wolframheizfadens 
unter verschiedenen Betriebsverhaltnissen 


Von 
BRANKO BERKES, Zagreb 


Mit 2 Textabbildungen 
(Eingegangen am 10. Juli 1960) 


Ubersicht: Wird der Wolframheizfaden in einer Réhre unter verschiedenen Betriebsverhaltnissen (kon- 
stante Heizspannung bzw. Heizstrom, konstante Fadentemperatur oder Anodenstrom) betrieben, so ent- 
stehen durch die verschiedene Verdampfungsgeschwindigkeiten Differenzen in der Lebensdauer des Heiz- 
fadens. Die Lebensdauer sowohl als auch die Differentialgleichung fiir den zeitlichen Verlauf der Halbmesser- 
verringerung und dadurch die Funktion ,(é) sind fiir alle vier Falle ermittelt worden. Die quantitativen Re- 
sultate kénnen durch den Gebrauch der in der Tabelle 1 angegebenen Daten benutzt werden. Fiir die Span- 
nungs- und Stromstabilisatoren sind auBerdem die Stabilitatsformeln (Anderung der stabilisierten Gré8e pro 
Stunde), welche fiir die praktische Anwendung von Sattigungsdioden in der Regeltechnik von Interesse sind, 
errechnet worden. 


Das Bediirfnis nach hochstabilen Strom- bzw. Spannungsdiskriminatoren fiir Stabilisations- 
gerate hat in dem letzten Jahrzehnten auch auf die Verwendung der Sattigungsdiode mit Wolfram- 
heizfaden gefiihrt [1], [2]. Abgesehen von geringer Empfindlichkeit auf Umgebungstemperatur- 
schwankungen und magnetische Streufelder und hohen Verstarkungsgrad zeigte sich die Satti- 
gungsdiode als besonders stabil im Zeitraum von mehreren Stunden und Tagen [2], [3]. Diese 
Eigenschaft kann jedenfalls nicht unbegrenzt ausgeniitzt werden, sondern wird von der Ver- 
dampfung des Wolfram stark beeinfluBt. Diese hangt aber von den Betriebsverhaltnissen des 
Heizfadens ab, d.h. ob der Wolframdraht bei konstanter Heizspannung, Heizstrom, Anoden- 
strom oder Temperatur verwendet wird. Jede von den erwahnten Betriebsarten verlangt ge- 
sonderte Betrachtungen und die quantitativen Resultate der Rechnungen werden in folgenden 
Abschnitten dargelegt. 

Wir bezeichnen M, in gcm™?s 7? als Verdampfungsgeschwindigkeit eines Einheitszylinders 
aus Wolfram, womit fiir die verdampfte Menge g in g eines geraden idealen! Fadens vom Halb- 
messer 7 in cm und Lange / in cm im Vakuum in kleinen Zeitintervall gilt 


dg =— M,:2r-1-dt. (1) 
Anderseits ist (y = 19,35 g cm 3 spez. Gewicht des Wolfram) 


dvi 27,2 1-7dr (2) 
und durch Vergleich mit (1) 
dy a _ My 
cae (3) 


Die Verdampfungsgeschwindigkeit kann nach [4] in kleinen Temperaturintervallen als 


(IT Temperatur des Heizfadens in Grad Kelvin) 


T \4 
Mi = Muo(7] 


ausgedriickt werden, worin sich die Indizes # auf einen beliebigen Anfangspunkt beziehen. Der 
Exponent yu ist von T,, abhangig, kann jedoch in Intervallen von je 100° K als konstant aufgefaBt 
werden [4], [6]. Damit nimmt die Differentialgleichung (3) die Gestalt 


Gee cage ie (zy (4) 


welche fiir alle Betriebsarten giiltig ist. 
Fiir konstante Temperatur des Fadens ist in jedem Zeitpunkt T = T,, und aus (4) folgt 


mit 7) als Anfangshalbmesser. 


1 Als ideal wird der Faden mit konstanter Temperaturverteilung iiber seine Lange bezeichnet. 
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Wenn die Lebensdauer der Diode wie friiher definiert wird, so kann man aus (8) fiir einen 
bestimmten 7), die endgiiltige Temperatur 7, errechnen (s. Tabelle). 

Um jetzt die Differentialgleichung fiir den zeitlichen Verlauf des Halbmessers zu gewinnen, 
mu 7/T,, in (4) durch (8) ersetzt werden, womit folgt 


3H 
dy a me (ay 


dt TN: 
und nach Auflésung (vy kann innerhalb von 100° K konstant gehalten werden) 


v—a 
y =(1 3u +—-2 t ee 
= : 


% 


(9) 


—¢« 20 Lr 


Die Lebensdauer des Fadens bei konstantem J F ist dann aus 


3u+ty—« 
| Bee ae y—o eee te 
Ty r 1 0,95 3u+yv—a ( 0) 


errechnet und tabelliert, nachdem alle Exponenten als Funktionen der Temperatur bekannt 
sind. Es ist ersichtlich, daB Ly, < L, ist, was physikalisch den Grund im — durch verkleinerten 
Halbmesser bzw. gro8eren Widerstand des Fadens — Zunehmen der Heizleistung hat. 
Ganz ahnlich wird die Endtemperatur und Lebensdauer fiir konstante Heizspannung er- 
rechnet, und wir erhalten aus 
Pile) C2F =) ei = au 
Up 


it TY 
Hig \ 


oder, nach Einsetzen in (4), den zeitlichen Verlauf der Halbmesserverringerung 
v+a 


fat gece Lane ta : gaa te (12) 


ae pra 20 Lr 


Die Lebensdauer bei U, = const ist aus (12) ermittelt zu 


bse 
Ly, = i * 20° ae Vict ie ) yta 


(13) 


(PSS 
In diesem Falle ist Ly, >L,, weil sich bei der Verringerung des Halbmessers auch die vom 
Faden verbrauchte Leistung vermindert. 
Das Verhalten von Wolframkathoden bei konstanter Heizspannung, Heizstrom oder Tempe- 
ratur ist einfach durch Messung des Anodenstromes nachzupriifen. Die Ergebnisse solcher 


Messungen mit konstanter Heizspan- 


nung wurden vor gewisser Zeit verffent- 
licht [9], [10], [12] und sie sind in guter 
Vereinbarung mit theoretischen Resul- 
taten aus dieser Arbeit. 

Die drei beschriebenen Betriebsarten 
— deren physikalischen Eigenschaften 
bisher schon bekannt, aber noch nicht 
mathematisch behandelt waren — sind ees 
vom Interesse, wenn die Speisung des ( Siittigungsdiode) 
Fadens von einer gesonderten AuBen- Bild 1. 
quelle hervorkommt. Soll aber eine 
Diode mit Wolframfaden in Stabilisationsgeraten als Diskriminatorelement eingebaut werden, so 
ist keiner von den dargelegten Fallen vorhanden. In einem Stromstabilisator nach Bild 1. [2] mit 


Stellglied 


x ° 0 Archiv fii 
62 B. BerxeS: Lebensdauer eines Wolframheizfadens — Eieke seat 


vorausgesetzter unendlich groBer Verstaérkung des Verstarketeiles ist in jedem Augenblick der 
Anodenstrom der Diode konstant. Aus der Emissionsgleichung [4], [6] folgt 


a ees a7 T. [%\* & bs tld on 4 
: aS, r= (") ¢ 2+ =| b= 52400 (14) 
und mit (4) gilt fiir kleine Temperaturanderungen 
F 
eA NE 
ae ( A 20 = ‘ (15) 


Die Lebensdauer findet man ahnlich wie in friiheren Fallen aus 


+d 
Day ee el EON ig (16) 


uta 

Fiir einen Stromstabilisator mu8 man wissen, in welchem Umfang sich der zu stabilisierende 
Strom durch Fadenverdampfung andert, d. h. mit welcher Stromstabilitat kann man bei einer 
gewissen Anfangstemperatur rechnen. Wiederholt vom Gleichgewicht der zugefiihrten und aus- 
gestrahlten Energie und der Formel fiir den Fadenwiderstand ausgehend folgen die Differential- 
gleichung 


dl dk d (ali dk (eli le d 
eer ba Scale eae CEA: f— 4 Ze (17) 
bzw. 
dl dT dy 
p=, 
2—! = (y— a) — = 


und mit (14) 
(18) 


welche Losung die Verbindung zwischen dem Heizstrom (stabilisierter Strom) und dem Halb- 


messer ergibt 
3A—v+o 


=(Z)~ . (19) 


% 


Es ist noch nachzupriifen, inwiefern man eine Verstarkerstufe von endlicher Verstéarkung K 
durch eine mit unendlich groBer Verstarkung ersetzen darf. Wegen einer durch dr verursachter 
Anderung dJ, entsteht im Verbraucherkreis eine Stromdifferenz 


di, =—dI, R,-§-K —_, 
Ry t 


mit € Koeffizient des ScHotrKy-Effektes [2], [7], &, Widerstand an der Ausgangseite des Regel- 
organs, S Steilheit des Regelorgans, die fiir S R, S> 1 mit Hilfe des RicHARDsoNschen Gesetzes in 


di - R O\ Gee dy b 
ool ee OK ee ee Noe pee oe 


iibergeht, und mit (17) bekommt man 


dT Ry BN; 
Pg) te 


Sry 


Sy Sy OF 
Ole Tate 


Das zweite Glied in beiden rechteckigen Klammern ist fiir die in der Praxis vorkommende 
Werte &, und J, stets vernachlassigbar und wir erhalten die Gl. (14) 


@ CO aie edie 
i ane i 
ip 


d. h. die Voraussetzung J, = const ist vollkommen gestattet. 
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Auf gleiche Art kann die zeitliche Stabilitat eines Spannungsstabilisator [1] (Bild 2) ermittelt 
werden. Es folgen die Gleichungen 


adUy _ 


Veron =| (20 
Che I a= (20) 
bzw. mit (14) 
2 UT OY ier ol Si wl A 
U; A Y 
und es 
sires Diskriminator 
Utp v 
unter gleichen Voraussetzungen wie im ao 


Falle des Stromstabilisators. Der zeit- 
liche Verlauf der Halbmesserverringerung und die Lebensdauer sind in beiden Fallen aus (15) 
und (16) zu bestimmen. 

Fiir Anwendung der Sattigungsdiode in Stabilisiergeraten ist statt der Lebensdauer mehr 
die Stabilitat? des Stromes bzw. der Spannung wissenswert. Wegen des hohen Verstarkungs- 
grades der Sattigungsdiode ist dieselbe allein fiir die Stabilitaét maBgebend und aus (19), (21) und 
(15) erhalt man annadhernd 


iE. —vy+ I a i) 
Ses SEE Ne an Slee al teats t=) (22) 
At 2A 20 Lr h 


fiir Stromstabilisatoren bzw. 


AU Atyv+e U ; V 
See ae — aie in (*) (23) 
fiir Spannungsstabilisatoren. Es ist selbstverstandlich, daB die durch (22) und (23) ausgedriickte 
Stabilitaét die 4uBerste Grenze darstellt, da die eigentliche Stabilitat wegen verschiedenen Ur- 
sachen [11], [12] statistischen Anderungen unterworfen ist. 

Der Verfasser spricht hiermit Herrn Dr. K. STEIMEL, Leiter der Forschung der Firma AEG, 
fiir freundliche Mitteilung tiber die Ergebnisse der Langzeitsstabilitat von Sattigungsdioden im 
Stromstabilisator und Kollegen Dipl. Ing. V. RADEKA und Dipl.-Ing. H. Basic fiir wertvolle 
Hinweise wahrend der Vorbereitung dieser Arbeit seinen aufrichtigsten Dank aus. 
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S. 304, Bild 2b und 
S. 305, Bild 4b: Der stehende Vektor hat jeweils den Wert 


V2 Jn @—y—o) 
o (1 + 0,) (1—@) 


der mit 2 w, rotierende Vektor den Wert 


t 


V2 Inve 74 
(1 +0,) (1 —9)" 


S. 304, GI. (17): Das Glied mit sin w, t in der eckigen Klammer muB ein m Mipnieeeienen haben. 


Ae 

<s 

S32. Gl,(61)* Innerhalb der dem sin y, als Faktoren in eckigen Klammern beigéordneten FS 
Ausdriicke miissen jeweils Summen und nicht Differenzen gebildet werden. — 
Ferner mu8 in dem Glied mit sin y, im Nenner o? erscheinen. as” 


= 


5: 383,-Gk (64): In dem in eckigen Klammern stehenden Ausdruck miissen die Glieder_ — 


mit Ky, und Ky, durch Pluszeichen verbunden werden. Das Glied mit he 
mu8 die Zeitkopstante 7, statt Ty, erhalten. 


S952. (68)): Hier muB in der ersten eckigen Klammer (1 + o) durch (1 — @) erset t 
werden, ferner im ersten Glied der zweiten Zeile o, durch o. Das letzt 


in runden Klammern stehende Glied dieser Formel ist durch ¢~?4/7@ au 
ersetzen. 


=. 314, 7, Zéile v. 0: 
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